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I.  INTRODUCTION 


A description  of  a standard  man  (the  Computer  Man) ^ has  been  built 
by  the  Target  Assessment  Branch  (TAB),  Vulnerability/Lethality  Division 
(VLD) , Ballistic  Research  Laboratory,  using  the  properly- chosen , hori- 
zontal, anatomical  cross  sections  (slices)  collected  by  Eycleshymer  and 
Schoemaker^  (see  Figure  1).  These  slices  have  been  uniformly  gridded  so 
that  the  physical  detail  of  the  Computer  Man  can  be  referenced  by  a 
Cartesian  lattice.  Current  descriptions  assume  that  all  tissue  in  a cell 
is  the  principal  tissue  in  that  cell.  The  two  currently-used  lattices 
have  a 0.5  cm  and  1.0  cm  horizontal  grid  respectively  and  both  descrip- 
tions use  2.6  cm  thick  lower-body  slices  and  1.2  cm  thick  upper-body 
slices. 

A phantom  of  the  man,  when  used  for  incapacitation  studies,  is 
assumed  to  be  standing  in  a box  that  barely  encloses  him.  The  box  is 
placed  in  the  first  octant  of  a Cartesian  coordinate  system  so  that  the 
sides  of  the  box  are  parallel  to  the  coordinate-axis  planes.  A bottom 
corner  of  the  box  is  located  at  the  coordinate-system  origin  so  that  the 
phantom  faces  the  positive-y  direction. 

A three-dimensional  matrix  is  constructed  whose  elements  reference 
the  Computer  Man  cells.  The  matrix  is  currently  filled  with  the  mortal- 
ity rankings  gathered  by  Cooper^  but  it  can  be  used  to  store  other, 
user-chosen  quantitative  descriptions.  The  indices  of  the  cell  enclos- 
ing a point  in  the  phantom  man  may  be  calculated  and  used  to  retrieve 
the  corresponding  matrix  element. 

The  FORTRAN  Computer  Code  RAYMAN  can  generate  rays  whose  origins  and 
directions  are  picked  from  user-selected  distribution  functions  and  pro- 
ject them  into  a phantom  of  the  Computer  Man.  Promenades  are  conducted 
along  the  rays  as  they  penetrate  tissue  and  a point  is  picked  uniformly 
along  each  step.  The  mortality  ranking,  M,  corresponding  to  the  cell  en- 
closing the  point,  is  retrieved  and  added  to  those  retrieved  for  the  same 
tissue  penetration,  R,  along  earlier  rays.  A histogramic  representation 
of  the  mean  mortality  ranking,  M,  versus  R is  calculated  at  the  end  of  a 
statistically  meaningful  sampling  of  rays  and  the  results  are  printed  in 
tabular  form.  This  output  may  be  used  to  assess  personnel  vulnerability 
to  the  ballistic  threat  as  discussed  by  Kokinakis  and  Bruchey.*^ 

^'Charles  A.  Stanley,  Midhael  Brown,  "The  Computer  Man,"  to  he  published 
as  a BEL  Report. 

2 

‘Albert  C.  Eyaleshymer,  Daniel  M.  Sohoemaker,  "A  Cross-Seation  Anatomy, " 
D.  Appleton-Century  Company,  Ina.,  New  York,  London  1938. 

3 

'Walter  R.  Cooper,  William  Kokinakis,  "Vulnerability  Rankings  of  Human 
Tissue, " to  he  published  as  a BRL  Report. 

4 

'William  Kokinakis,  William  J.  Bruchey,  "An  Engineering  Approach  to  the 
Assessment  of  Personnel  Vulnerability ,"  American  Defense  Preparedness 
Association,  October  1975. 
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II.  CODE  CAPABILITY 


Simulated  projectiles  (rays)  can  be  started  from  origins  whose 
locations  are  picked  from  a user-selected  choice  of  nine  different 
distributions.  The  distributions  are: 

1.  User-supplied  coordinates.  The  present  dimensions  of  the 
storage  arrays  permit  the  user  to  submit  a maximum  of  100  points. 

2.  Points  generated  uniformly  on  the  overhead  hemisphere.  The 
user  will  choose  the  radius  for  all  hemispherical  origins. 

3.  Points  generated  uniformly  on  the  underfoot  hemisphere. 

4.  Points  generated  uniformly  on  the  X > 0 hemisphere. 

5.  Points  generated  uniformly  on  the  X < 0 hemisphere. 

6.  Points  generated  uniformly  on  the  Y > 0 hemisphere. 

7.  Points  generated  uniformly  on  the  Y < 0 hemisphere. 

8.  Points  generated  uniformly  on  a horizontal  circle  parallel 
to  the  plane  on  which  the  man  is  standing.  The  user  will  provide  the 
height  and  radius  of  the  circle. 

9.  Points  generated  on  an  overhead,  quarter  circle  arc.  The 
plane  of  the  circle  will  include  the  Z-axis  and  make  a user-supplied 
aspect  angle  with  the  XZ  plane.  Points  will  be  picked  from  a user- 
chosen  distribution.  The  two  possible  distribution  functions  are,  (a) 
uniform,  and  (b)  the  sine  of  the  polar  angle  of  the  radius  vector  of 
the  point  (see  Figure  2). 

The  directions  taken  by  the  rays  are  picked  from  a user-selected 
choice  of  three  different  angular  distributions.  These  distributions 
are : 

1.  Isotropic  about  the  point  of  origin.  In  practice,  the 
direction  of  the  ray  is  usually  chosen  so  that  it  falls  within  the  solid 
angle  subtended  by  the  sphere  enclosing  the  man  box.  A normalization 
factor,  that  is  the  quotient  of  the  solid  angle  subtended  by  the  sphere 
divided  by  4''^,  is  calculated  and  printed  in  the  output. 

2.  Bivariate  normal  about  the  line  from  the  ray  origin  to  a user- 
supplied  point  within  the  man.  This  capability  was  provided  to  the  code  to 
simulate  the  marksmanship  of  an  individual  aiming  a weapon  toward  a point 
on  the  man.  The  horizontal  and  vertical  standard  deviations  of  the 


marksman's  efforts  in  degrees  are  provided  by  the  user. 


3.  User-supplied  direction  cosines.  The  code  will  automatic- 
ally restrict  the  problem  to  one  history  per  point  when  this  option  is 
exercised. 

The  size  of  the  array  required  to  accommodate  a matrix  of  a com- 
plete description  exceeds  the  memory  size  of  many  computers.  This 
limitation  is  circumvented  by  dividing  the  man  into  horizontally-bounded 
sections  whose  descriptions  are  of  a manageable  size  and  storing  the 
matrix  of  each  section  as  a record  in  a disc  file.  The  matrix  of  a 
section  will  be  chosen  and  loaded  into  memory  as  needed  by  RAYMAN  gener- 
ated commands  and  promenades  will  then  be  conducted  through  that  part  of 
the  phantom.  RAYMAN  stores  the  necessary  tracking  parameters  when  a 
promenade  along  a shotline  enters  an  inactive  section  (a  section  whose 
matrix  is  not  currently  loaded  in  the  array  MAN) . The  parameters  of 
such  a latent  promenade  will  be  retrieved  at  a later  time  for  further 
tracking.  This  technique  reduces  the  number  of  time-consuming  disk 
loading  operations  as  compared  to  the  number  required  when  tracking  each 
shotline  to  its  termination  before  generating  another  ray. 

The  logic  flow  involved  in  generating  a ray,  conductine  a nromcnadc 
along  a shotline,  and  retrieving  and  tallying  scores  is  now  discussed. 

More  detail  can  be  obtained  from  the  flow  charts  of  the  main  routine 
(Figure  4)  and  principle  subroutines  (Figures  5-9).  A complete  listing 
of  the  program  is  given  in  the  Appendix.  A description  of  the  subroutines 
and  the  identification  of  the  principle  variables  is  given  in  Section  III. 

A ray's  origin  and  direction  are  picked  from  the  chosen  density 
functions.  The  ray  is  projected  toward  the  phantom  man  and  subroutine 
BOX  is  called  to  determine  if  the  computer  man  box  is  intercepted.  The 
length  SLR  of  the  segment  intercepted  is  calculated  by  the  subroutine 
for  those  rays  intercepting  the  box.  Subroutine  TRACK  is  then  entered 
and  a promenade  is  started  along  the  shotline  at  the  entrance  point  of 
the  ray  if  the  section  of  the  man  being  entered  is  the  section  whose 
matrix  is  currently  stored  in  the  array  MAN.  The  promenade  step  width 
is  set  to  the  histogramic  bin  width  when  a promenade  is  started.  The 
parameters  of  the  promenade  are  stored  in  latent  storage  if  the  matrix 
stored  in  MAN  does  not  belong  to  the  section  being  entered. 

The  promenade  is  continued  inside  the  box  until  tissue  is  encountered. 
The  distance  traversed  through  the  box  PB0X  is  incremented  after  each  step 
and  compared  to  SLR  to  determine  if  another  full-width  step  in  the  box  is 
possible.  The  distance  to  a section  crossing  SLRC  is  calculated  by  sub- 
routine CR0SS  to  determine  if  another  full  width  step  in  the  active 
section  is  possible  (see  Figure  .5). 

The  first  score  is  made  when  the  promenade  enters  tissue.  The  score 
is  estimated  stochastically  when  the  complete  length  of  that  step  is  in 
tissue.  This  is  acconq^lished  by  picking  a point  with  equal  probability 
at  any  site  along  the  step  and  retrieving  the  M-value  from  the  array  MA>.' 


I 

SLR  = [(X2-X,)2  * ( Z2-Z,)^]2 


P{x»yiZ)  = Current  Position  Along  Promenade. 
PBDY  = PBDY,  ♦ PBDYj 


o « Scoring  Point. 


* Sampling  Step  Width. 


Figure  3.  A Cross  Section  of  the  Computer  Man 
that  Contains  a Shotline 
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Figure  7.  Latis  Subroutine 
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Figure  8.  Trace  Subroutine 
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Figure  9.  Track  Subroutine 


fjr  the  cell  enclosing  that  point.  This  score  is  then  added  to  the  first 
bin  of  the  array  VI.  The  square  of  the  score  is  added  to  the  first  bin 
of  the  V12  array.  The  array  NVI,  used  to  tally  the  number  of  promenades 
penetrating  that  amount  of  tissue,  has  its  first  bin  incremented  by  one. 
The  amount  of  tissue  penetrated  during  the  current  promenade  PBDY  is  set 
to  the  histogramic  bin  width  DELR. 

The  score  is  calculated  exactly  when  the  promenade  initially  enters 
tissue  and  the  entering  step  is  only  partially  in  tissue.  The  final 
score  deposited  in  the  Vl-array  will  also  be  made  for  a tissue  penetra- 
tion depth  of  DELR  but  more  than  one  promenade  step  will  be  needed  to 
attain  that  penetration.  A partial  score  SCPT  is  obtained  for  the 
tissue-entering  step  by  calculating  the  product  of  the  fraction  of  a bin 
width  intercepted  by  a cell  and  the  M-value  of  that  cell  and  summing 
those  products  over  all  cells  encountered  along  the  step. 


This  is  accomplished  in  the  code  by  calling  subroutine  TRACE  to  cal- 
culate the  necessary  cell  intercepts  and  to  retrieve  the  necessary  M- 
values.  The  intercept  lengths  are  stored  in  DIST(l,m)  and  the  M-values 
in  DIST(2,m).  A partial  score  SCPT  is  then  calculated  that  is  given  by 

NCELL 

SCPT  = ^ ^ DIST(l,m)  • DIST(2,m)  (1) 

m=l 

where 

NCELL  = the  number  of  cells  intercepted  during 
the  step 

m = the  running  index  over  cells  intercepted. 

The  variable  PBDY  will  be  incremented  by 

NCELL 

PBDY  = 

where 

H[DIST(2,m)  - e]  = 


The  variable  PB0X  will. 

Subroutine  EVEN  is 

the  additional  step(s)  needed  to  obtain  a full  bin  width  of  tissue  pene- 
tration. These  latter  score(s)  may  be  estimated  stochastically  or  cal- 
culated analytically  as  determined  by  the  tissue  or  non-tissue  status  of 


5"^  DIST(l,m)  • H[DIST(2,m)  -e] 
m=l 


(2) 


the  Heaviside  step  function, 
a positive  infinitesmal . 

as  usual,  be  incremented  by  the  full  step  width, 
now  called  to  calculate  the  partial  score(s)  of 


i 

f 


i 
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successive  cells.  The  total  score  for  a full  bin-width  of  tissue  pene- 
tration is  added  to  the  VI  array  and  the  other  variables  and  arrays  are 
properly  incremented.  The  step  width  is  then  reset  to  DELR  and  a normal- 
gaited  promenade  through  tissue  is  ready  to  be  started. 

That  pace,  conducted  in  subroutine  TRACK,  is  maintained  until  the 
promenade  exits  the  box  or  current  section,  or  leaves  tissue.  A point 
is  picked  uniformly  along  each  step  and  the  M-value  associated  with  the 
cell  enclosing  that  point  is  retrieved  from  the  array  MAN  and  stored  in 
the  proper  bin  of  the  Vl-array. 

Subroutine  TRACE  is  called  again  when  checks  predict  that  the  next 
step  will  exit  tissue.  The  variable  values  retrieved  or  calculated  by 
TRACE  are  used  to  calculate  a partial  score  SCPT,  as  described  in  Eq.  2 
for  the  amount  of  tissue  penetrated  during  the  step.  Subroutine  EVEN 
is  then  entered  and  the  promenade  is  continued  until  tissue  is  re- 
entered (or  the  box  is  exited  or  a new  section  is  entered) . The  partial 
score  for  a segment  of  a length  sufficient  to  complete  the  bin  is  cal- 
culated when  tissue  is  reencountered  and  is  added  to  the  previously  cal- 
culated partial  score.  The  total  score  is  added  to  the  proper  VI  bin, 
all  variables  are  updated,  and  control  of  the  promenade  is  returned  to 
subroutine  TRACK  where  a normal  gait  is  resumed. 

The  partial  score  will  always  be  calculated  for  the  last  step  of  a 
promenade  in  the  box  (regardless  of  whether  that  step  was  encountered 
in  subroutine  TRACK  or  EVEN).  The  partial  score  is  calculated  in  the 
customary  manner,  normalized  to  a full  bin  width,  in  a statistically 
valid  way,  and  added  to  the  proper  VI  bin. 

The  partial  score  will  always  be  calculated  for  the  last  step  before 
a section  crossing.  SCPT  is  calculated  in  the  usual  manner  and  placed  in 
latent  storage  along  with  the  other  necessary  parameter  values.  Sub- 
routine LATENT  performs  this  storage  as  well  as  the  later  retrieval. 

The  parameter  values  for  100  promenades  can  be  stored  as  latents.  j 

The  flag  ILATX  is  used  to  direct  the  storage  and  retrieval  of  these  para-  j 

meters.  ILATX  is  set  to  1 at  the  beginning  of  each  editing  unit  of  j 

histories.  It  remains  at  this  value  until  latent  storage  is  saturated  | 

or  all  promenades  required  by  the  current  editing  group  have  been  started.  : 

This  value  of  1 directs  the  code  to  store  latents  when  necessary  but  to  i 

not  yet  retrieve  latents.  j 

ILATX  is  set  to  2 when  the  parameters  for  the  100  latent  promenades  | 

have  been  stored  and  will  remain  at  that  value  until  all  histories  re-  j 

quested  by  the  current  editing  group  have  been  started.  This  mode  orders 
the  retrieval  of  the  next  latent  promenade  that  is  entering  the  currently 
loaded  section  of  the  man  when  the  tracking  of  the  previous  promenade 
(either  latent  or  virgin)  ended  in  a section  crossing,  with  a subsequent 
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latent  storage.  A box  exit  of  the  previous  promenade  will,  on  the  other 
hand,  return  a space  for  the  potential  storage  of  a new  latent  that  is 
derived  from  a presently  unstarted  promenade.  New  rays  may  therefore  be 
generated  and  their  promenades  conducted  until  one  of  them  has  to  be 
stored  as  a latent. 

The  sweep  of  the  latent  storage  arrays  for  retrieving  latent  prom- 
enades entering  the  currently  loaded  section  of  the  man  is  started  when 
ILATX  is  set  to  2.  The  Vl-matrix  of  the  section  having  the  most  latent 
promenades  making  an  entrance,  is  loaded  into  MAN  at  the  completion  of  a 
sweep  and  a new  sweep  will  be  initiated. 

ILATX  is  set  to  3 when  all  histories  in  an  editing  group  have  been 
started.  New  rays  will  not  be  generated  and  latent  promenades  will  be 
retrieved  and  tracked  until  latent  storage  is  exhausted.  Restorage  of  a 
latent  promenade  after  its  retrieval  and  tracking  is  possible  in  this  as 
well  as  the  preceding  mode  of  ILATX. 

The  editing  of  results  may  be  implemented  in  the  ends  of  each  his- 
tory, at  the  end  of  each  group  of  histories  emanating  from  a point,  or 
at  the  end  of  a run.  The  mean  of  the  score  and  its  standard  deviation 
are  calculated  for  each  bin.  These  values,  as  well  as  the  number  of 
promenades  penetrating  that  amount  of  tissue,  are  printed  in  tabular 
form  along  with  identifying  captions. 


III.  SUBROUTINES  AND  VARIABLES 


Subroutine  Descriptions 


1.  Driver  Routine.  This  routine  reads  in  the  user-supplied  data 
and  instructions.  Once  calculated  values  of  variables  are  derived,  be- 
fore entering  the  ray-generating  cycle,  the  ILATX-routing  to  different 
modes  of  operation  and  choosing  of  Vl-matrices  for  loading  into  the  array 
MAN  are  also  conducted  in  this  routine. 


2.  Subroutine  AD.  This  routine  directs  the  logic  flow  to  the 
user-chosen  ray  angular  distribution. 


3.  Subroutine  ARC.  This  subroutine  picks  an  origin  for  the  ray  on 
an  overhead,  quarter-circle  arc.  These  points  may  be  picked  from  two 
distributions.  The  distributions  are:  (a)  uniform,  and  (b)  sine  of  the 
polar  angle  of  the  point. 


4.  Subroutine  BIVAR.  This  subroutine  picks  the  direction  of  a 
ray  from  a bivariate  normal  distribution  about  the  line  from  the  ray 
origin  to  a user-supplied  point  usually  (but  not  necessarily)  located  in 
the  man.  The  vertical  and  horizontal  standard  deviations  of  the  marks- 
man are  supplied  by  the  user. 
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5.  Subroutine  BOX.  The  man  is  assumed  to  be  barely  enclosed  by 
a rectangular  box.  This  subroutine  will  project  the  current  ray  toward 
the  box  and  calculate  the  length  SLR  of  the  segment  intercepted.  The 
coordinates  of  the  entering  point  of  the  ray  are  also  returned  to  be 
used  as  the  start  of  the  promenade. 


6 .  Subroutine  CHECK, 
valued  on  the  cell  boundaries. 


Tlie  Computer  Man  description  is  double- 
Subroutine  CHECK  inspects  rays  to 


determine  if  they  will  travel  on  these  boundaries.  A ray  that  would  have 
such  a path  will  have  the  critical  coordinate  of  its  origin  randomly 
translated  by  a small  amount  so  that  its  path  will  not  lie  on  a cell 
boundary . 

7.  Subroutine  CR0SS.  This  subroutine  is  used  to  calculate  the 
distance  from  the  beginning  of  the  current  promenade  step  to  the  nearest 


z boundary.  This  distance  SLRC  is  returned  to  the  main  program. 


8.  Subroutine  EDIT.  This  subroutine  directs  the  calculation  of 
the  mean  and  standard  deviations  of  accumulated  scores.  It  will  also 
print  the  results  in  tabular  form  along  with  identifying  captions. 

9.  Subroutine  EVEN.  This  subroutine  is  called  by  subroutine 
TRACK  after  a promenade  has  entered  (or  reentered)  tissue  or  has  been 
revived  after  storage  as  a latent.  The  promenade  gait  will  be  brought 
into  phase  with  the  histogramic  bin  structure  of  tissue  penetration. 
Subroutines  SCORE  and  TALLY  are  called  in  this  subroutine. 

10.  Subroutine  HEMI . This  subroutine  picks  a ray-origin  point 
uniformly  on  any  of  the  user-chosen  hemispheres. 

11.  Subroutine  INKS.  This  subroutine  calculates  the  distance 
from  a point  on  a shotline  in  a cell  to  the  exit  of  the  shotline  from 
the  cell.  This  distance,  as  well  as  the  exit  point,  are  returned  to 
subroutine  TRACE  for  use  in  calculating  partial  scores. 

12.  Subroutine  ISQT.  This  subroutine  is  used  to  pick  the  direction 
of  a ray  uniformly  about  its  origin.  In  practice,  the  wasteful  practice 
of  picking  directions  for  rays  that  would  not  intercept  the  man  is 
improved  by  picking  directions  that  will  project  a ray  within  the  solid 
angle  subtended  by  a sphere  barely  enclosing  the  man  box.  A normalization 
factor  that  is  the  ratio  of  this  solid  angle  to  Att  is  calculated  so  that 
the  user  may  calculate  absolute  probabilities.  The  direction  parameters 
are  converted  to  direction  cosines  and  the  latter  quantities  are  returned 
to  the  main  program. 

13.  Subroutine  LATENT.  This  subroutine  is  called  to  store  or 
retrieve  latent  promenades.  The  0 or  1 value  of  the  first  variable 

IL  in  the  argument  directs  toward  a retrieval  or  a storage  respectively. 
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14.  Subroutine  LATIS.  This  subroutine  locates  the  cell  that  en- 
closes the  point  PI.  The  x-axis  index  (II),  y-axis  index  (Jl),  and 
z-axis  index  (Kl)  are  returned.  The  I,  J,  and  K indices  are  derived 
from  II,  Jl,  and  Kl  and  are  used  in  their  place  for  retrieving  scores 
so  that  program  abortion  caused  by  truncation  approximation  will  be 
avoided  at  the  box  boundaries.  A point  on  a shot  line  that  lies  on  a 
cell  boundary  will  be  assumed  to  lie  within  the  cell  being  entered. 

15.  Subroutine  LOAD.  The  principal  function  of  this  subroutine 

is  to  load  the  Vl-matrix  of  a chosen  section  into  the  array  MAN  from  disc. 
A secondary  role  is  to  call  subroutine  TEST  to  generate  a test  MAN  array 
when  requested  by  the  user. 

16.  Subroutine  LOCATE.  The  role  of  this  subroutine  is  to  determine 
the  section  of  the  Computer  Man  that  includes  the  z-axis  index  KP 
(absolute) . This  index  is  then  converted  to  its  value  relative  to  the 
bottom  of  that  section  and  returned  as  KP  (relative) . 

17.  Subroutine  SCORE.  Tliis  subroutine  retrieves  the  score  from 
MAN(I,J,K).  The  score  is  returned  in  both  fixed  point  and  floating  point 
mode.  It  is  anticipated  that  this  subroutine  may  be  rewritten  as 
required  by  other  fortlicoming  descriptions  of  the  Computer  Man. 

18.  Subroutine  SORC.  This  subroutine  provides  the  routing  to  the 
chosen  ray-generating  routines.  It  will  also  calculate  the  distance 
from  the  ray  origin  to  the  center  of  the  box  or  to  the  aiming  point  in 
the  man . 

19.  Subroutine  TALLY.  This  subroutine  will  add  the  score  that 
was  estimated  or  calculated  to  that  already  accumulated  in  the  proper 
bin  in  the  VI  array.  The  bin  of  the  array  that  tallies  the  number  of 
shotlines  reaching  that  penetration  in  tissue  NVI  is  also  incremented 
by  1.  The  square  of  the  score  is  added  to  the  VI2  array  to  be  used  at 
the  completion  of  the  editing  group  to  calculate  standard  deviations. 

20.  Subroutine  TEST.  This  subroutine  is  used  to  build  a test 
pattern  of  scores  in  the  MAN'  array.  It  can  be  easily  rewritten  to 
satisfy  the  needs  of  the  user. 

21.  Subroutine  TRACE.  This  subroutine  is  used  to  advance  a 
distance  R along  a shotline  while  calculating  the  distances  penetrated 
in  intercepted  cells  and  retrieving  the  scores  in  these  cells.  The 
number  of  cells  in  which  the  segmeiit  R lies  is  also  calculated  and 
returned  as  NCELL. 

22.  Subroutine  TRACK.  This  subroutine  initiates  the  promenade 
along  a newly-generated  shot  line  or  resumes  the  promenade  along  a 
retrieved  latent  shotline.  Subroutine  EVEN  is  called  \i;hen  needed  to 
bring  the  promenade  in  phase  with  the  histogramic  bin  structure.  TRACK 
will  conduct  an  even-gaited  pace  at  other  times  (except  for  box-,  section-. 
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or  tissue-exits)  with  stochastic  estimation  of  scores.  TRACK  will  also 
calculate  partial  scores  for  the  noted  exceptions  to  the  even  gait  when 
the  preceding  step  was  in  TRACK. 

23.  Subroutine  TRANS.  ITie  Computer  Man  box  is  assumed  to  be 
placed  entirely  in  the  first  Cartesian  octant  with  a corner  located  at 
the  axis-origin.  The  RAYMAN  subroutines  that  pick  ray  origins  will  use 
quarter-circles  or  spheres  that  are  centered  at  the  origin.  This  sub- 
routine will  translate  the  ray-origins,  after  being  picked,  so  that  they 
will  lie  on  circles  of  spheres  centered  at  the  center  of  the  man  box. 

24.  Subroutine  VAR.  This  subroutine  calculates  the  standard 
deviation  of  each  scoring  bin. 

B.  Variable  Description 

1.  BULSI ( 3) . The  coordinates  of  the  aiming  point  in  the  Computer 
Man.  It  is  used  only  when  the  bivariate  normal  angular  distribution  is 
used. 


2.  BKPT.  The  height  at  which  the  thickness  of  the  slices  changes. 

3.  DELR.  This  is  the  width  of  the  histogramic  bins  in  which 
scores  are  accumulated.  DELR  is  also  the  preferred  step  width  of  the 
promenade  through  tissue. 

4A.  DIST(2, IZ) . The  score  retrieved  in  the  IZth  cell  tracked  by 
subroutine  TRACE. 

4B.  DIST(1 , IZ) . The  length  of  the  segment  intercepted  by  the  IZth 
cell  in  the  tracking  by  subroutine  TRACE. 

5A.  FACT(l)  1=1,10.  An  adjustment  factor  for  quantizing  a ranking 
in  the  Vl-matrix.  All  elements  in  F'ACT  are  currently  set  to  unity. 

5B.  FACT (11) . This  clement  of  FACT  is  set  to  1.0  and  is  used  in 
calculating  partial  scores. 

6.  GCH.  The  height  of  the  great  circle  above  the  plane  = 0. 

7.  GCR.  The  radius  of  the  great  circle  about  the  vertical  axis 
of  the  Computer  Man. 


8A.  GRID(I)  1=1,2.  The  x-  and  y-cell  width  of  the  Computer  Man. 
8B.  GRID(3) . The  z-cell  width  of  the  lower  part  of  the  Computer 

Man. 

8C.  GRID(4) . The  z-cell  width  of  the  upper  part  of  the  Computer 

Man. 


9.  The  x-axis  cell  index  of  the  Computer  Man.  This  quantity 

is  adjusted  at  the  box  boundaries  to  prevent  computer  truncation  approxi 
mations  from  causing  the  logic  errors  that  would  lead  to  incorrect 
answers  or  run  abortion. 

10.  I^.  ITie  same  as  1 for  internal  cells  but  it  is  not  adjusted 
at  the  box  boundaries. 

11.  IBKPT.  The  z-axis  of  the  top  slice  of  cells  possessing  a 
thickness  of  GR1D(3) . 

12.  IDIST(IZ) . The  fixed  point  score  in  the  IZth  cell  tracked  in 
subroutine  TRACE. 


13.  IHIST.  The  running  index  for  the  rays  generated  from  a point. 

14.  IHIT.  The  number  of  interceptions  made  in  the  walls  of  the 

Computer  Man  box  by  the  ray. 

15.  ILAT.  The  running  index  for  the  retrieval  of  latent  promenades. 

16.  1 LATX . The  flag  for  directing  latent  promenade  storage  and 

retrieval.  The  value  ILATX  = 1 directs  toward  storage  only.  The  value 

ILATX  = 2 permits  either  storage  or  retrieval.  The  value  ILATX  = 3 
permits  retrieval  only. 

17.  IPT.  The  running  index  for  the  ray  origin  point. 

18.  ISECT.  The  current  section  of  the  Computer  Man  stored  in  tlie 
array  MAN . 

19.  ISECTT.  Ihe  section  of  the  Computer  Man  in  which  a designated 
point  is  located. 

20.  ISEED.  A fixed  point  number  used  as  an  argument  b>'  the  (0,1) 
random  number  generator . 

21.  ISL.  The  maximum  number  of  horizontal  slices  of  the  Computer 
Man  stored  in  a section.  This  number  is  the  same  for  all  sections 
except  possibly  the  topmost  section. 

22.  ISC.  The  fixed  point  value  of  the  score  in  a cell. 

23.  J.  The  y-axis  cell  index  of  the  Computer  Man.  This  quantit\- 
is  adjusted  at  the  box  boundaries  to  allow  for  computer  truncation 
approximations . 

24.  The  same  as  J for  internal  cells  but  it  is  not  adjusted 
at  the  box  boundaries. 


25A.  JLAT(1 , IP) . The  section  number  into  which  the  IPth  latent 
promenade  is  entering. 


25B.  JLAT(2 , IP) . The  number  of  histogramic  bins  into  which  scores 
have  been  tallied  for  the  IPth  latent  promenade. 

26.  JSEEb . The  fixed  point  argument  for  generating  a pair  of 
random  numbers  picked  from  a normal  distribution. 

27.  JEXIT.  A flag  showing  that  a cell  of  subroutine  EVEN  has 
led  to  a box  exit  by  the  current  promenade. 

28.  1^  The  z-axis  cell  index  of  the  Computer  Man.  This  quantity 
is  adjusted  <it  the  box  boundaries  to  allow  for  computer  truncation 
approximations  during  calculations  along  the  promenade. 

29.  The  same  as  K for  internal  cells  but  it  is  not  adjusted 
at  the  box  boundaries . 

30.  MAN (I ,J ,K) . The  array  used  to  store  the  incapacitation  input 
data  for  the  current  section  of  the  Computer  Man. 

31.  NARC . A flag  for  choosing  the  distribution  from  which  an 
origin  point  of  a ray  on  the  quarter-circle  is  picked.  The  value  NARC  = 

0 leads  to  picking  by  a distribution  given  by  the  sine  of  the  polar  angle 
while  the  value  NARC  = 1 leads  to  picking  a point  uniformly  on  the  arc. 

32.  NCELL . Tlie  number  of  cells  penetrated  by  the  shotline  segment 
when  subroutine  TRACE  is  used. 

33.  NDEL . The  number  of  histogramic  bins  into  which  scores  have 
been  added  for  the  current  promenade . 

34.  NFACT.  A flag  for  determining  if  FACT  values  are  to  be  read 
in  as  part  of  the  input  or  if  the  values  are  to  be  defaulted  to  1.0. 

35.  NFLGl . A flag  for  determining  the  distribution  from  which 
ray  origins  will  be  picked.  The  options  are  discussed  in  the  user- 
supplied  input  section  of  this  report. 

36.  NFLG2 . A flag  for  determining  the  distribution  from  which 
ray  directions  are  chosen.  The  options  are  discussed  in  the  user- 
supplied  input  section  of  this  report. 

37.  NFLG3 . The  flag  for  determining  the  editing  group  to  be 
used.  The  options  are  discussed  in  the  user-supplied  input  section  of 
this  report. 

38.  NUIST.  The  number  of  rays  to  be  started  at  a point. 


L. 
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39.  NHITS . The  number  of  rays  intercepting  the  Computer  Man  for 
an  editing  group. 

40.  NLAT( ISECT) . The  number  of  latent  promenades  that  are  entering 
section  ISECT. 

41.  NLATT.  The  number  of  latent  promenades  stored. 

42.  N LOAD . Tlie  number  of  times  that  a different  Vl-inatrix  was 
loaded  into  MAN. 

43.  NPTS . The  number  of  points  from  which  rays  are  to  be  started. 

44.  NSECT.  The  number  of  sections  into  which  the  Computer  Man  is 
divided . 

45.  NTEST . A flag  to  determine  if  test  values  are  to  be  generated 
for  the  array  MAN. 

46.  NVI  (I)  1 = 1 ,200 . The  counter  for  determining  tlie  number  of 
shotlines  that  have  penetrated  a distance  in  tissue  greater  than  (1-1) 
•DEDT. 

47.  NX.  The  number  of  Computer  Man  cells  in  the  x-direction. 

48.  NY . The  number  of  Computer  Man  cells  in  the  y-direction. 

49.  The  number  of  Computer  Man  cells  in  the  z-direction. 

50.  P(I)  1=1 ,5.  The  coordinates  c>f  the  ray  origin. 

51.  PBDY ■ The  distance  in  tissue  penetrated  during  a promenade. 

52.  PB0X . The  distance  in  the  man  box  penetrated  during  a promenade. 

53.  PC(I)  1=1,3.  The  coordinates  of  the  point  at  which  a promenade 
will  cross  the  next  section  boundary. 

54.\  PLAT( 1 , IP) . The  initial  value  of  PBDY  for  the  IPth  Intent. 

54B.  PLAT(2 , IP) . The  initial  value  of  PB0X  for  the  IPth  latent. 

54C.  PLAT(5,IP) . The  value  of  SLR  for  the  IPth  latent. 

55.  PS  (I)  1 = 1,3.  The  coordinates  of  the  beginning  of  tiie  step  in 
a promenade. 

56.  PR(I)  1=1,3.  The  coordinates  of  the  uniformly-picked  point  on 
a promenade  step. 
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57.  PE (I)  1=1,3.  The  coordinates  of  the  end  of  the  step  in  a 
promenade . 

58.  PTAR(I,IP)  1=1,3.  The  coordinates  of  the  beginning  of  the  IPth 
latent  promenade  or  the  coordinates  of  the  origin  of  a user-provided  ray. 

59.  PTAR(I,IP)  1=4,6.  The  direction  cosines  of  the  shotline  or 
ray  described  in  item  58. 

60.  PI  (I)  , P2(l)  , 1 = 1 ,5.  The  coordinates  of  the  intersection  of 
the  ray  with  the  X=0  and  X=XL3  plane  respectively.  This  usage  occurs 
in  subroutine  BOX. 

61.  P5(I) , P4(I)  , 1 = 1 ,3.  The  coordinates  of  the  intersection  of 
the  ray  with  Y=0  and  Y=XL2  plane  respectively.  This  usage  occurs  in 
subroutine  BOX. 

62.  P5 (I)  , P6 (1)  1 = 1,3.  The  coordinates  of  the  intersection  of 
the  ray  with  the  2=0  and  Z=XL3  plane  respectively.  This  usage  occurs 
in  subroutine  BOX. 

63.  PI  (1)  , P2( 1)  . These  arrays  are  also  used  in  subroutines 
TRACE  and  INX3  to  store  the  coordinates  of  the  start  or  entrance  of  a 
segment  into  a cell  and  its  termination  or  exit  from  the  cell. 

64.  PDEL . Tlie  amount  of  tissue  penetrated  in  a new  bin.  This 
variable  is  used  in  subroutine  TRACK. 

65.  PLl . This  variable  is  used  in  subroutine  TRACE  and  INX3. 

The  plane  X=PL1  includes  one  side  of  the  cell  which  could  be  the  exit 
boundary  of  the  shotline  from  the  cell. 

66.  PL2 . This  variable  is  used  in  subroutine  TRACK  and  TRACE. 

The  plane  Y=PL2  includes  another  side  of  the  cell  which  could  be  the 

exit  boundary  of  the  shotlinc  from  the  cell. 

67.  PL3 . This  variable  is  used  in  subroutine  TRACK  and  TRACE. 

The  plane  Z=PL3  includes  a third  side  of  the  cell  which  could  be  the 

exit  boundary  of  the  shotline  from  the  cell. 

68.  I^.  The  radius  of  the  sphere  that  barely  encloses  the 
Computer  Man  box. 

69.  The  distance  from  the  ray  origin  to  the  center  of  the  man 
box.  It  is  also  used  as  the  distance  from  the  ray  origin  to  the  point 
BULSI . 

70.  SCPT.  A partial  score  that  is  calculated  exactly  when  a 
promenade  is  being  brought  into  phase  with  the  histogramic  bin  structure. 


71.  SDl . The  standard  deviation  (degrees)  in  the  vertical  direction 
of  the  shots  triggered  by  a marksman. 

72.  SD2 . The  standard  deviation  (degrees)  in  the  horizontal  dir- 
ection of  the  shots  triggered  by  a marksman. 

73.  SLR.  The  length  of  the  segment  of  the  ray  intercepted  by  the 
Computer  Man  box. 

74.  SPCT.  The  ratio  of  the  solid  angle  subtended  by  the  Computer 
Man  box  sphere  to  4it. 

75.  SLRC.  The  distance  from  the  beginning  of  a step  to  a section 
crossing.  This  variable  is  defaulted  to  2,000,000  when  W(3)  = 0.0. 

76.  SPLT.  An  array  for  storing  the  heights  of  the  section  boundaries. 
The  values  are  stored  in  ascending  order. 

77.  The  floating-point  score. 

78.  VI (IP) . The  bin  for  accumulating  VI  scores  whose  penetration 
in  tissue  is  between  (IP-1)*  DELR  and  IP*  DELR. 

79.  VI2 (IP) . The  array  for  storing  the  square  of  the  scores 
deposited  in  the  VI  array. 

80.  W(IZ) , IZ=1 ,5 . An  array  for  storing  the  direction  cosines  of 
the  current  ray . 

81.  WP (IZ) , IZ=1 , 5 . An  array  for  storing  the  direction  cosines  of 
the  line  from  the  ray  origin  to  either  the  riming  point  (BULSI)  or  the 
center  of  the  Computer  Man  box. 

82.  XLl . The  plane  X=XL1  is  an  outside  boundary  of  the  Computer 
Man  box. 

83.  XL2 . The  plane  Y=XL2  is  an  outside  boundary  of  the  Computer 
Man  box . 

84.  XL3 . The  plane  Z=XL3  is  an  outside  boundary  of  the  Computer 
Man  box . 

85.  Z9 . The  length  of  the  part  of  a segment  of  the  shotline  that 
is  in  the  current  cell.  This  variable  name  is  used  in  subroutine  TRACE 
and  INX3. 
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IV.  RAYMAN  INPUT 


1.  TITLE 

A inaximum  of  80  alphanumeric  identifying  characters. 

(8AI0) . 

2.  ISEEU  , JSEEU 

ISEED  = (0.1)  random  number  routing  argument. 

JSEED  = normal  random  number  routine  argument  (2110). 

3.  NTEST,  NFACT,  NARC 

NTEST  = Flag  for  generating  test  numbers  for  the  MAN  array. 

NTEST  = 1,  generates  numbers  for  MAN. 

NTEST  = 0,  loads  data  from  disc  to  MAN. 

NF'ACT  = Flag  for  filling  the  FACT  array. 

NFACT  = 0,  reads  in  user- chosen  numbers. 

NFACT  = 1,  all  elements  in  FACT  are  defaulted  to  unity. 

NARC  = Flag  to  determine  distribution  of  ray  origins  on  the  quarter 
circle. 

NARC  = 0,  points  are  picked  from  the  sine  of  the  polar  angle, 
distribution. 

NARC  = 1,  points  are  picked  uniformly  along  the  arc. 

(1015) 

4.  NFLGl,  NFLG2,  NFLG3,  NPTS,  NHIST 

N'FLGl  = Ra>’  origin  flag. 

NFLGl  = 1,  read  in  points. 

= 2,  uniformly  distributed  on  the  upper  hemisphere. 

= 3,  uniforml)'  distributed  on  the  lower  hemisphere. 

= 4,  uniformly  distributed  on  the  front  hemisphere. 

= 5,  uniformly  distributed  on  the  back  hemisphere. 

= 6,  uniformly  distributed  on  the  right  hemisphere. 

= 7,  uniformly  distributed  on  the  left  hemisphere. 

= 8,  uniformly  distributed  on  the  great  circle. 

= 9,  distributed  on  quarter-circle  arc. 

NFLG2  = Ray  direction  flag. 

NFLG2  = 1,  isotropic  angular  distribution. 

= 2,  bivariate  normal  angular  distribution. 

= 3,  read  in  directions. 


28 


NFLG3  = Editing  group  flag. 

NFLG3  = 1,  composite  editing  group. 

= 2,  each  point  editing  group. 

= 3,  each  history  editing  group. 

NPTS  = Number  of  ray  origin  points . 

NHIST  = Number  of  rays  per  point . 

(1015) 

5.  NX,  NY,  NZ,  IBKPT,  ISECT,  NSECT 

NX  = Number  of  x-axis  cell  boundaries  in  box. 

NY  = Number  of  y-axis  cell  boundaries  in  box. 

NZ  = Niimber  of  z-axis  cell  boundaries  in  box. 

IBKPT  = Number  of  z-axis  cells  up  to  slice  thickness  change. 

ISECT  = Initial  section  of  the  man  loaded  into  MAN. 

NSECT  = Number  of  sections . 

(1015) 

6.  ISL 

ISL  = Number  of  slices  per  section  (except  for  possibly  the  topmost 
section) . 

(1015) 

7.  GRID(I),  1=1,4  DELR 

GRID(l)  = x-axis  cell  width  (cm)  . 

GRID(2)  = y-axis  cell  width  (cm) . 

GRID(3)  = z-axis  lower  cell  width  (cm) . 

GRID(4)  = z-axis  upper  cell  width  (cm)  . 

DELR  = promenade  step  width  (cm) . 

(6E10.3) 

8A.  IF  NFLGl  = 1,  provide 

P0(I),  1=1,3;  WC(I),  1=1,3  IF  NPTS=1 

PTAR(I,IPT)  1 = 1,6  IF  NPTSM 

P0(1)  or  PTAR(1,IPT)  = x-coordinate  of  ray  origin 

P0(2)  or  PTAR(2,IPT)  = y-coordinate  of  ray  origin 

P0(3)  or  PTAR(3,IPT)  = z-coordinate  of  ray  origin 

W(l)  or  PTAR(4,IPT)  = x-axis  direction  cosine  of  ray 

W(2)  or  PTAR(5,IPT)  = y-axis  direction  cosine  of  ray 
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W(3)  or  PTAR(6,IPT)  = z-axis  direction  cosine  of  ray 
(6E10.3) 


8B.  IF  NFLGl  = 2-7,  provide 
R2 

R2  = radius  of  hemisphere 
(6E10.3) 

8C.  IF  NFLGl  = 8,  provide 
GRC,  GQl 

GRC  = Radius  of  circle  (cm) . 

GQl  = Height  of  circle  above  z=0 

(6E10.3) 

8D.  IF  NFLGl  = 9,  provide 
R2,  THETA 

R2  = Radius  of  arc  (cm) 

THETA  = azimuth  angle  (degrees)  with  x-axis 
(6E10.3) 

9.  IF  NFLG2  = 2.  provide 
SDl,  SD2 

SDl  = azimuthal  standard  deviation  (degrees)  of  the  marksman. 
SD2  = elevation  standard  deviation  (degrees)  of  the  marksman. 

(6E10.3) 

10.  IF  NFLG2  = 2,  provide 
BULSI(I)  1=1,3 

BULSI(l)  = x-coordinate  of  aiming  point 
BULS1(2)  = y-coordinate  of  aiming  point 
BULSI(3)  = z-coordinate  of  aiming  point 

(6E10.3) 

11.  IF  NFACT  = 0,  provide 
FACT(I)  1=1,10 

FACT(l)  = adjustment  for  VI  = 1 

FACT(IO)  = adjustment  factor  for  VI  = 10 
FACT(ll)  = 1.0  always 

(6E10.3) 
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V.  TEST  PROBLEM 

The  input  needed  for  a test  problem  is  given  in  Table  I.  This  data 
is  presented  in  the  format  in  which  it  is  printed  in  the  output  by  the 
computer.  The  output  average  Vl-Scores  of  this  test  problem  is  presented 
in  Table  11. 


VI.  CONCLUSIONS 

The  RAYMAN  Computer  Code  has  been  run  successfully  on  the  BRLESC 
computer.  Readers  desiring  a copy  of  the  code  may  send  a blank  magnetic 
tape  to  the  Target  Assessment  Branch,  Vulnerability/Lethality  Division, 
ATTN:  DRDAR-BVL,  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland  21005.  Users  are  cautioned  that  the  BRLESC  names  of  certain 
library  routines  may  not  agree  with  those  used  by  other  computer  facilities 
and  will  need  to  be  changed. 


! 

i 

1 
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Table  I.  Sample  Input 


onooF  00  v-r.H  fct  • O.o’OOE  oo 


iNCRE.'lkNT 


Table  II.  Sample  Output 

/IJL  Tioex  OFV  'IQ  hays 


1* 

<1 

9 

0,70'^'OC 

nu- 

0.3iiOE 

01* 

9 

* 

20 

* 

9 

0,7150E 

■lU- 

0 , 3 3 4 2 E 

01  * 

9 

* 

20 

3* 

9 

0.7I50E 

.)U- 

■3.2  75^E 

)I* 

9 

* 

20 

<** 

* 

9 

^,73C0e 

U*- 

'J.2  5I3E 

11* 

9 

* 

20 

5* 

9 

r,72^lE 

3 I 

.J.2-.5''E 

u* 

9 

* 

20 

6* 

* 

9 

0,6  76SP 

:U- 

0 , 20  30E 

01* 

9 

* 

17 

7* 

9 

, 63  F dE 

ou- 

0,34  l3E^ 

01* 

9 

* 

17 

8* 

:ai 

9 

0,6647E 

JU- 

C,3635E 

u* 

9 

* 

17 

9-4 

4I 

9 

^,72C0E 

01  + - 

■J.  27  3 IE 

•Jl* 

9 

♦ 

15 

10* 

« 

9 

'3.63576 

OU- 

0.3?  )'^E 

•)1* 

9 

* 

14 

II* 

* 

9 

■'>,P077E 

' # i ♦ * 

G,22)C'E 

•'a* 

9 

* 

13 

12* 

9 

0.3231E 

JU- 

0,2242E 

01* 

9 

* 

1? 

13* 

<E 

9 

'',315hE 

1 ♦ - 

0.244^E 

•n* 

9 

* 

13 

I't* 

9 

3,34l7E 

cu- 

i'.2'72E 

a* 

9 

* 

12 

15* 

9 

0,73iaE 

0 1-f- 

U.2'72E 

a* 

9 

* 

1 1 

16» 

« 

9 

0,7571E 

3-t*- 

0,3452c 

01=, 

9 

* 

11 

17* 

9 

0.7323E 

.u+- 

0,3‘'2‘^E 

a* 

9 

* 

10 

IS* 

« 

9 

3.  ni  IE 

Oi.*- 

1.325  7E 

01* 

9 

* 

C. 

19* 

<! 

9 

'•',77C3E 

•u>- 

0. 3''3o£ 

0 1 * 

9 

* 

c 

20* 

9 

O.aorOE 

0 1 * 

0 , 1 « 2 7 E 

n* 

9 

« 

21* 

♦ 

9 

0,0250E 

: 1 f- 

G.  173JE 

a * 

9 

« 

r 

22* 

9 

0,3V-5  0E 

I ■♦•  - 

0,2  n='E 

■01^ 

9 

* 

r 

23* 

« 

9 

3,'3322E 

ii+- 

,3I47e 

a* 

9 

* 

k 

2<V* 

•At 

9 

0,3571P 

.u- 

j,  113^E 

01* 

9 

« 

7 

25* 

# 

9 

n,^2«liF. 

n+- 

0. 1 OOE 

a* 

9 

« 

7 

2t)* 

9 

9 

^,•*7146 

0I>- 

0,407  IE 

a» 

9 

f 

27* 

9 

9 

''.*357E 

i:  I +■  - 

■J,  30  )eE 

a* 

9 

7 

2 e * 

9 

9 

>.7714E 

01*- 

l,3  3-J2E 

n* 

9 

* 

t 

29* 

9 

9 

2, 7571  E 

/ i +- 

G , 3 ' 39  c 

01* 

9 

« 

7 

30* 

9 

9 

0,7:jC0E 

J.363 'E 

a* 

9 

« 

31* 

9 

9 

. 30('0E 

i * - 

0.2  7 3'^E 

a* 

9 

u 

•32* 

9 

9 

”3,  '>0''OE 

lU- 

u . 30  P 2 r 

a* 

9 

* 

z 

33* 

9 

9 

'•>,7u<'0c 

>u+- 

0. 3’''  Jbe 

a* 

9 

* 

2 

lA* 

9 

9 

6 7r9E 

n*- 

0 , 4 ’ '7  1 E 

a* 

9 

♦ 

■5 

35* 

9 

9 

■3.  lO'-OE 

52  + - 

O.O’'  luE 

lr..i 

9 

* 

36* 

9 

9 

3,  lOt'OE 

173  >^£-06* 

9 

* 

2 

37* 

9 

9 

^.lOCOE 

32*- 

0.9-^30E 

a * 

9 

* 

* 

3t)* 

9 

9 

■3.  lUf'OE 

12  + - 

ii . 9 0 3 G E 

Ol* 

9 

♦ 

1 

39* 

9 

9 

0,  lOroE 

Jl  + - 

').9000E 

a* 

9 

* 

1 

3 J* 

9 

9 

O,2O''0E 

■,l*- 

0.00  90 E 

Ol* 

9 

* 

1 

AI* 

9 

9 

0.2O'‘05 

11*- 

v,9o9CE 

a* 

9 

« 

1 

'♦2* 

9 

9 

O.SjroE 

U*- 

0,909('E 

a* 

9 

* 

1 

33* 

« 

9 

''.AjroE 

•:  l+- 

0 , QTOUP 

a* 

9 

* 

1, 

'»4* 

9 

9 

^.50''0E 

•'  1*- 

0 , 90  OOP 

a * 

9 

* 

1 

35* 

9 

9 

^ ,60r3E 

31*- 

J.99  3CP 

31  * 

9 

XI 

1 
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DMEmSION  TITLe(8» 

COi1UN/nAH>J/M*M|  110, 55*6), NTEST 

CniMnN/L*TEN/ JLAT ( 2, l0Ul,Nt*T(20),IL*T,NL*Tr, I SECT, lL/kTX,ISECTT, 

• PL*T(4,J00|,PTAR(6, 1 0 0 1 , PH nT , PHOX , 8CP T , NOEL , I H I ST , I PT , NARC , MPT S 
COHHiJN/LATI»/OHin<A),RXPT,ISXPT,SPLT(lO| 
C0HHliN/P3l>iT/PS)3),PR(JI,PE(3),pC(3),DELR,‘)Fl.03 
C0>1HUN/H|s:/Rl,R2,3CR,aCH,NFL3l,NFLa2,SDl,SU2,Qwl,8H3 
C0>iHjr4/r^AC/SLR,l(l.t,XL2,XL3,NHIT3,DIST(2,10)  , IDIST)  10),NCELL 
COHHJ>)/E VAL/V I (20U),Nvl(20a),VI2(2Cn),FACr(lll 
COiHUN/S3RS/Pp(3),H(3)>P(3),POI3),BUl.SI (3), SPOT- ISEED,JSEED, 

• SnT,CST 

C0'iK;-j/i.AT/HX,NY,NZ,NSeCT,I9L 

COiMus/LDOE/RUOAO 

500  FORrt4T(Fl0.1) 

501  FORHaTIO >101 

502  FORHATdOISI 

503  FORHAri6El0.3) 

504  FORHATdH  ) 

505  FORMAT ( 1 Ml ) 

506  F0RMAT(2d0l 

507  F0RMAT16( AlO, |5,5X) ) 

508  FORMA  r ( 4<  4 1 0,  Ed  .4,9X1  ) 

509  F0RMAr)2A10,A3,ed,7) 

510  FORM4T(A10,A4, 15) 

511  F0RMAr(A10,ll(F9.4,lX)) 

512  forma  r(2IA10,d0,2X)) 

513  FORMAT! 19F10.5) 

514  FORMATd9X,12,4X,F10.5) 

515  format ( 6(E20. 14 ),/, 120) 

516  format ( 31 F20. 14)1 

REAoi3,5on  title 

REAU(b,506)  ISFeD,JS£EU 

READI5,502)NTEST,NFACT,NARC 

READ  Co, 502)  NFLSI,NFL32,nFL83,NPTS,RH13T 

IF! HFL32.F3.3)  NHIST  ■ 1 

READ  I 5, 50?)  NX,RT,NZ,IBKPT,ISECT,NSECT 

READI5,5U2)  I 3L 

RE  AO ! 5, 503  ) OR  1 0, UELR 

IF ! NFLSt .ST, 1 ) 00  TO  6 

IF ! NPTS.  IT. 1 ) 00  TO  4 

READ! 5,503)  PO, 4 

OOTO  9 

6 IF1NFL31.lt. 9)  30  TO  9 
HEAD!5,5Q3)  R2, theta 
theta  ■ 3.  141S<>»THETA/180.0 
CST  ■ COSFITHETA) 

SNT  • SlNFlTHCTA) 

00  TO  9 
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BBr  AVVJUCIE 


4 4CA0(9,503I  ( (PT4R(I > JI > ("t >6) « Jat >NPTS> 

SOTO  9 

9 IFINFLSl.OT.TI  90  TO  S 
9EaD(9.909I  92 
80  TO  9 

B 9EAn(3,909)  OCRiOCH 
>|RITE(6>S04) 

9 IF(NFU32.NE.2I  00  TO  JO 

RE*D(9,9a3l  SU1,S02,0>«1>0H2 

RE4D(9,903)  RULSl 
10  lF(NFACr.E9.0|  OOTO  99B 
DO  997  (■1,10 
FACT (II  ■ 1.0 
997  CONTINUE  S OOTO  996 
990  CONTINUE 

read(5,903)  Fact 
996  CONTINUE 

FACTItll  ■ 1.0 

9KPT  • OR] P( 3)»FL0ATF ( IBKPT) 

XL1«QRID( 1 l•FLOATF(NX) 

XL2"09IU(2l«F|.0ATF(Nr) 
XL3«o'(RT*8Rm(4|»FL0ATF(NZ-IBKI'TI 
R1  • (0.5»XLll»»2  ♦ (0.9«XL2I»»2  ♦ ( 0 . 9 • Xl.3  I • *2 
R1  ■ 3URTF(R1| 

J>0 

DO  777  l»l,N9ECT 
J-J+ISl, 

|F( J.OI.IBXPT)  00  TO  776 
XJ»FL0ATF( J) 

SRLT( I l-XJ^OKIDIO) 
liO  TO  777 

776  XJ«FL0ATF( J-le<PTI 

SRLT(l|a9KFT»(XJI*GRIDI4l 
777  continue 
ILATX  • 1 
•«RITEI6,909) 

DO  770  I-I.NSECT 
uRI TE(5, 51  A I ( I .SRWTI  I I I 
778  CONTINUE 

»RITE(6,50|I  TITLE 
nRITEIB.SOAI 

CAPl  ■ lONISEEO  ■ iCAP2  ■ lOHJSEEO  • 
i4R|Tkl5.9t2)  CAPI,ISEEU,CAP2.JSEE0 
WRITE  I 6, 504  I 

CAPl  ■ 10MNTE9T  ■ 1CAP2  ■ lOMNFACT  • 

•iR|T£I6.507ICAP1.NTEST>CAP2,NFACT 

>RITE(6,504I 

CAP!  • IOHNARC* 

nRITEI  6. 507  I CAPl .NARC 

WRITEI 6.5041 
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Btsi  AVAllAB'i  COPY 


c*i>l  • IQH  SQ  PTS  > IC*P2  ■ IOMNO  hist  • 

■<RlTtt»,S07IC*Pl,SPTS,tAP2,NH|ST 

« R I T E ( 6 , 3 0 4 I 

C*»l  • 10hIU*TX  • 

hr  1 Tt 1 6, so; ) CAP  I , I LATX 

hRI Tt ( i,  504  I 

CAPl  • IOmnELGI  • JC‘P2  ■ 10HNEL52  • ICAPJ  • IOhNFLQS  ■ 
«RlTt(b.S0;)CAPl,^FLGl»CAp^,HFL3^,CAP3,NFL03 

*R I TE ( 6,  304  I 

CAPl  • lOHNG  X'S  • JC*P2  • IOHNO  T'S  ■ ICAPS  ■ lOHNO  Z'S  ■ 

HRITtlb,50;)C4Pl,NA,CAR2,NY,CAP3,NZ 

"RITE 6*5041 

CAPl  • 10H18XPT  • IC»P2  • lOHSECT  NO  ■ SCAP3  ■ lOHNO  SECT  • 
CAPA  • lOHNO  SLICE  ■ 

hR1TE(6,50;>Capi,1bKPT.CAP2,ISECT,CAP3.NSECT,CAP4,ISL 

HR  I TE I S,  504  I 

CAPl  • lOHX-CEUL  ■ $CAP2  • lOHY-CELL  ■ SCAP3  ■ IOmZI-CELL  • 
CAPA  ■ lO-tZZ-CELU  • 

HRIT£(5*508)CAPl,QR[Din*CAP2,0RlD(2l,CAP3,GHID(3),CAP4,SRlD(4» 

HR  1 TE I 6*  504  I 

CAPl  H lOHR-iNCREH  • 

hR ITE ( 6*  SOSICAPI , OELH 

hRI TE I 6*  504  1 

CAPl  • lOHSOURCE-X  "ICARZ  • lOMSOURCE-Y  hICAPS  ■ IOHSOURCE-Z  ■ 
hRITE |b*50SICAPl,P0ll|*CAp2,PUI2),CAP3*P0l3) 
nR I TE I 6*  504  ) 

CAPl  H lOHA-DIR  CS  •SCAP2  • IOHY-DIR  CS  "SCAPS  b IOMZ-UIR  C8  « 

«RITE(5*50#)CAP1,h(1),CaP2*h(2)*0AP3*H(3I 

hRI TE I 6*  504  1 

CAPl  B 10H30X.RAD  ■ SCAP2  b 10HPT«RAD  b 

hRITE(5*505ICAPI,R1*CAP2*R2 

*)R  I TE  1 6*  504  I 

CAPl  B lOHX-SlDE  B ICAP2  ■ lOHY-SIDE  • SCAP3  b IOHZ-SJDE  b 

HRirE16*509)CAPl,XLl*C*P2*XL2,CAP3,XL3 

HR  1 TE ( 6*  504  ) 

»RI TE I 6*  504  ) 

CAPl  • lOHFACT  HAT  B 
bR I TE ( 6*  51 1 ) CAPl , fact 
hRITE(6*504) 

CAPl  B lOHiELT  RAD  b*C*P2  ■ lOMBELT  HOT  b 

hRITE(5*50S)CAPI*OCR*C*P2,OCH 

hR I TE ( 6 * 30  A ) 

CAP)  B lOHSULSI(l)  BSCAP2  a 10hPJLSI(2)  bjcaps  a lOHBULSKS)  a 
BR|Tt(6*505)CAPl,«OLSlll),CAP2*BJLSI(2l*CAP3*BULSl(3) 

HR  I Tt ( 6*  504  ) 

CAP1b10m3R£AT  CIRC*  CAi'2»lOHLE  SN-*NpL*  CAP3b3HE» 


>S  CAPAb 

hR I TE I 6 
hRI TE  I 5 
bR 1 TE ( 6 


IOHLE  CS-ANGL 
502)CAP1*CAP2,CAP3,SNT 

504  ) 

509)  CAP1*CAPA,CAP3*CST 


.y9 


f 


c 

KOlL 


7 

1 

1 1 
t2 

13 

14 

3 

20 


[ 

t 

994 


21 

995 


102 

103 


loa 

113 
1 1 4 


»ll»ITE(5,S04) 

C4P1  ■ lOHHOU  SD  • *C*f*2  • lOMVERT  SO  ■ 

■«9ITE(5,508)C*P1,SU1,C*P2,SD2 

49 1 TE ( 8, 504  I 

C4P1  ■ 10MH-S4  rCT  •IC4P2  • 10MV-9P  FCT  ■ 
4B|TEt8,508)  C4PlfQ»‘l,CAP2»0H2 
4P I Tb ( 8, 504  I 
4P1TE(8,504) 

*RITb(8,504l 
«P I TE I 6, 504  I 
«RIT£(8,504I 
t«RITE(  8,504l 

C*U1.  LOAOIISECTI  $ NUOAD-l 

RLATT  • 0 

IPT  ■ 0 

IHIST  • 0 

IPT  • IPT  ♦ I 

IHI3T  4 ImiST  ♦ I 

CAUL  30RC 

CAUL  dOXI IdITl 

IPI IHIT.UT.2)  80  TO  995 

CAUL  CHECK 

IP(NPL33.LT.3)  30  TO  21 
CAUL  LATISIP, I,J,KI 
CAUL  LOCATEIK,  I3ECTTI 
CAUL  LATENT! 1 , IPT, P) 
nLATT«NUATT*1 

nLAT(ISECTT)4NLAT(ISECTT)41 
ISECT  ■ JLATI 1 , IPT) 

IF( ISECT. Ea.O)  GO  TO  99S 

CAUL  L0A0(I3ECT)  S NLQAD  ■ NLOAO  * 1 

CAUL  TRACK(O) 

30  TO  994 
call  TRACKll) 

IF ( NPL33.GT.  1 ) SO  TO  103 
IF( IhIST.LT.NhISTI  SO  10  108 
IF( IPT.Ea.NPTSl  GO  TO  102 
MIST  4 0 3 IPT  4 IPT  ♦ 1 

SO  TU  lOS 

ILATX  43  3 GO  TO  113 

I F 1 NFL33.GT .2  I GO  TO  118 
1F( MIST.LT.NHIST)  go  10  108 
30  TO  102 

IF(NLATT.LT.100I  go  TO  14 

IF( ILATX. E3, 2)  GO  TO  109 

ILATX  4 2 

ILAT  4 0 

ILAT  4 iuaT  ♦ 1 

IF( JLATI 1 , ILAT) .NE. ISECT)  00  TO  120 
call  TRACK(O) 
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IF( ILATX.FB.3)  GO  TO  120 
109  IF( JL*T ( I , ILAT 1 .EU.O I 00  TO  14 
120  IFl ILAT.UT.NLATT)  00  TO  114 

115  NUAT( 13ECTI  ■ 0 
l-l*  lbECT»l 

2 I*I>1 

IF(NLAT(I).3T.NLAT(ISECTin3ECT«I 

IF(  1 .LT. NSECT lOOTO  2 

IF ( NLAT ( 1 SECT ) .EU.O ) GU  TO  lift 

CAUL  LOao(ISECT)  1 NUUAD  ■ NLOAO  ♦ I 

30  TO  113 

116  CAPl  - lOH'tO  HAN  HITS  CAP2  ■ 4HS  ■ 
kPITE(6>310)CAp1,CAp2,NV|(1) 

"9  1 r£( 6>  504 ) 

CAPl  ■ 1 OHNOHHAL I 2ATJ  CAP2  • lOHION  FACT0RtCAP3  ■ 3H  ■ 

■RITE(5,509ICAP1,CAP2, LAPS, SPOT 

NR1Tc(6,504I 

CAPl  ■ 10H13EED  • »CAP2  ■ lOHjsEEO  • 

"R|Te(6,512»  CAPl, 1 SEEK, CAP2,J SEED 
"R 1 TE ( 6, 504  I 

CAPlalOHNO  DISC  LOS  CAP2"4HA0  ■ 

"R I TE ( 6, 51 0 I CAPl ,CAP2, NLOAD 
CAUL  Eutr 
ILATX  ■ I 

NLATT  "01  NHITS  ■ 0 
220  IF( NFU33.E3. 1 I 00  TO  117 
IF(NFUi3.0T.2l  GO  TO  999 
IF( IPT.EJ.NPTS)  30  TO  117 
GO  TO  12 

999  IF( IhIST.LT.nhIST)  00  10  M 
IF( IPT.LT.NPTSI  go  to  12 

117  STOP 

98765  CONTINOE 

98766  CONTINUE 

"RITE (6, 51  7 I SCM|,SCATI,nlOaD,NTE3T,IL*T,NLATT, I SEC T , I L A TX , I SECTT 

517  FORXAT  ( 2( £12. 4, 2X I , 7 ( 1 4, 2X ) » 

"»ITE(6,518I  PHD»,PrtOX,3CPT,NO£L,IHlST,IPT,NARC,NPTS 

518  FONH* T ( 3 ( El  2. 4, 2* • , 51  ! 6, 2X n 
•*RITt(6,519)  ORIO.HXPT/IHKPr.SPLT 

519  FO"M*ri5(E10.2,2X),I4,t20.12,2X,/,6(E20.12,lX),/,6(E20,12,lX),/, 
*6  1 E20. 1 2,  I X ) I 

«RITE(5,520)  PS,pR,PE,PC,DELR,NFL03 

520  F0RH*T19tE12.A,lX),/,4(ei2,4,lX),I4l 

"RITt(6,521)  R1,R2,GCR,OCH,NFLG1,NFL32,SD1,SB2,OH1,9H2 

521  F0yH4TI4(El2.4,lX),2(l4,2xl,4(E12,4,Ul) 

WR|TE(6,522I  SLR,XL1,XL2,XL3,NmITS,01ST, I DI ST.NCELL 

522  FOR-AT  (4(E12.4,ix),I5»/,1U(E12.*,1X),/,10(E12.4,1X),/,11(I5,1X)( 
»R1TE(6,523)  FACT 

523  FORMAT ( 1 1 ( E12.4  I ) 

PRI TE ( 5, 524  I pP, ", P,PO, SULS  I ,SPCT,I3EED,JSEED 
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S24  F0i«H<T(6(ei2.4,|X),/>6IEt2.4,tXI,/,4(C12.4,U)>2l  I9.IXI> 


i<RlTe(&,  504l 
wB|TE(5,525l  NL*T 
525  FOKMxr ( 10(2X» 1511 
«<9ITEIb,  504) 

«i«ITE<  5,525)  JLAT 
mRITE(6,S04I 
STOP 
END 
C 
C 

SUSROUriNE  *u 

CONHON/n{3C/Rl,R2,8CR,aCH,NFt.01>NPL02>SDl>SD2,6Ht,SN2 
|F(NFL82,3r.l)  OOTO  2 
CAUL  ISOT 
30  TO  ] 

2 CONTINUE 

IF(NFLS2.3T.2»  OOTO  3 
CAUL  SIVAR 

3 RETURN 
END 

C 

C 

SUBROUTINE  ARCtNARC) 

CON  NARC  a 1 IS  UNIFORM,  NARC  ■ 0 IS  ISOTROPIC  IE  COSINE  mEISHTED 
C0NN0N/S3RS/mP(3I ,w(3l ,P(3I ,P0(3),BUUSI (3),3RCT,ISEED,JSEED, 

• SnT,CST 

C0NMUN/NISC/Rl,R2,0CRiaCH,NFUQt ,NFUS2,SDt ,SD2,Gw1,3n2 

if(narc.eq.o)  go  to  2 

1 XlaUKANOH  ISEED) 

X2aUKAN31 ( ISEED) 

XIS  ■ Xl**2  S X2S  ■ X2**2 
T«X1S*X2S 

if(t.3t,i.o)  goto  » 

SNPHla2.0«Xl»X2/T 
CSPHI«A8SF((X1S-X23I/T1 
GOTO  3 

2 continue 

CSRHl  ■ URAN3l(ISEEDI 

SNPHI  • SORTFd.O  - C3PHI*«2I 

3 CONTINUE 
P(3)«R2»CSPHI 
Pin  aRB^SRPHI  aCST 
P(2)>R2*SNPHI*SNT 
RETURN 

END 

C 

C 

subroutine  BIVAR 

C0NN0N/S0RS/MPI3) tN(3l,PI3l .P0(3I (BUUSKSI ,SPCT,ISEED,JS£ED, 
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Btsi  availAd’^  copy 


• S«<T,CST 

COHMOS/h|sc/N|,H2,9CR>BCH,NFLQI,>IFLS3>SD1>SO3,QN1,0n2 

IFii«P13).Ea.»l.O)  GO  TO  9 

IF(mF< JI.EO.l.ni  GOTO  B 

IFt  Jt'KIi.Ea.U  GOTO  I 

PSI""H(3I 

PSI"*'<CCa8(P5l  I 

THET»«p(2)/80RTP(1.0-WP(3)'»«2» 

THET«ARCC03(TmETI 

IPRIn^J 

SPCT-l.O 

SDl  ■ 501*3. 14159/lSO.O 
SD2  ■ 502*3.14159/180.0 
1 call  f^^Avll  (Rnn1,RNN2,  JSEEOI 
0EL*SU2*RNR2 
RSN3-04AN31 t ISEEO) 

IF |RN 13.3T.0.5)UEU>*0EU 

PSlM*PSl*nEL 

N(3I*C33F(P8|H) 

0EL-30t*RNRl 
RNN3*JRA<13i  ( ISEEDI 
1F(«M  '(3.3T.0.5)DEL«-DEL 
TH£TM*rHET*DEL 
SI9T«534TF(1.0-WP<3)***) 

4(2l*3|Nr*C08F(THETH» 

« ( 1 1*51  NT *51 NF( THETHI 
4R1TEI 5, 10)  4P,W,P,RNNI,RNN2»TMETH 
10  FORPAT|6(E10.3>2XI  ) 

5 return 

ENO 

C 

c 

SUSROUTlRE  BOX(IHIT) 

DIRENSIOR  pi  9(3) ,P0UT( 3) ,NPL (21 

DMENSIOR  P1(3I»P2(3),P3(3),P4(3)»P5(3I,P6(3I 

CORNUN/S3R3/>iP(3l,N(3)rPI3)>PO<3I.BULSI ( 3) .SPOT. 1 SEED« JSEED. 

* SrT.CST 

C0RMON/TRAC/SlR»XLI,XL2,XU3,>JH1TS,D1STI2,10I,IDI8T(10|»NCEUL 

1HIT*0 

CON  PI  And  P2  are  X«0  AND  X»XLl  INTERCEPTS 
Pt I 1 1*0.0  I P2( 1 )*XL1 
IFIhI 1 I .EQ.O.QI  GOTO  1 
29*-P| 1 )/Rl  11 
P1(2I**(2)*29*P(2I 
PI I3I*4I3)*Z9*PI3I 
Z9*Z9*XLl/4( 1 1 
P2(2)*«i(2I*Z9*P(2» 

P2( 3) *4 ( 3) *Z9*P (3) 

CON  X*0  Ano  X*XLI  planes 

CON  CHECK  FOR  BOX  INTERCEPT 
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-IFIPI  (21 .SE.XLZ.OR.PI  (3) .LE.0.0I9OTO  9 
|P(P1  ( 31. OE.Xt.3. OP, PHJI.UE, 0.0)90X0  9 
|H|T>1  S NPLd)"! 

00  lOU  Ial,3 
PMt  I )»Pl  ( I ) 
too  COMTINOE 

9 |P(P2(2).GE.XL2. OR, P2(2).LE. 0.0)90X0  201 
IP(P2(3),0E.Xl3,0H,P2(3),|.E,0.O)3OXO  20 1 
IHIXaMIXat 
DO  1 Ial.3 
|P( MIX. EO, 2)00X0  6 
PM(l)aP2(l)  9 NPL(l)a2  t OOXO  7 
6 P0UX( 1 (■P2( I I $ NPL(2)a2 
1 CORXMUE 

201  in  MIX. EQ. 2)00X0  20 
COH  P3  XNO  P«  XRE  Yao  AND  Xa  XL2  MXCRCEPT9 
1 P3(2)a0.0  S P4(2)axL2 
IF(al2) .EO.O.O)  90X0  3 
iga.PI  2 ) /•((  2) 

P3( 1 )aa( 1 ) tZOaPI 1 I 

P3( 3)aa( 3) aZOaPI 31 
Z9aZyaXL2/R(2l 
P4( 1 )aa( 1 )*Z9»PI 1 I 
P4I  3 )aa I 3 ) aZ9«PI 3) 

C0*1  xao  AND  raXL2  PLANES 

IP|P3( 1 ) .OE.XUI .0R.P9( 1 ) .lE.0,0)  90X0  10 
IP(P3I3).OE.XL3,ON.P3I3),LE,0.0>  00X0  10 
Ml  XaMI  Xal 
00  lUl  |at,3 
IFI  MIX. E9. 2)00X0  9 
PMII)aP3(I)  t NPL(l)a3  S 00X0  101 
aOJX(  I )a(>3(  I ) s NPL(2)a3 
CONX  M JE 

IFI  P4( 1)  .SE.XLI .0R.P4( 1 ) .LE.O.O)  90X0  202 
|F|P4(3).gE,XL3,0R,P4|3).LE.0,0)  00X0  202 
MlXaMIXat 
00  12  Ml. 3 
IFI Ml X.EO, 2)00X0  II 
PMIMaPAll)  S NPL(l)a4  S 00X0  12 
PODXI I )aP4( I ) S NPU(2)a4 
CONX  M JE 

|F|  Ml X.EO. 2)30X0  20 

P9  AND  P6  ARE  Zao  AND  EaXL3  MXERCEPXS 
P9l3)a0.a  t P6l3)aXL3 
IFMID.EO.O.O)  30  XO  203 
iga.PI3)/N(3) 

PSI 1 )a4 I 1 ) aZOpPI  1 ) 

P9l2)a4|2)«Z9«P(2) 

Z9aZ9AXL3/4(3) 

P6I 1 )a4( 1 ) •Z9*P|  1 I 


9 

101 

10 


1 1 
12 
202 
CON 

3 
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1 


COM 


14 

15 


16 


IS 

203 


COM 


19 

50 


20 


21 


22 


23 

25 


51 

52 


M6(2)««( 2)*Z9«P(2> 

ZaO  ANU  £■(13  PLANES 

IP ( P5( 1 ) .OE.XLl .0N.P5( I ) .LE.0,0)  80T0  16 

ir(P5(  21 ,3E.XL2.0«.P5r£> .LE.O.OI  OOTO  16 

|HIT«Mir»l 

00  15  1-1,3 

1PI |H| r.E0.2)00T0  14 

PiNlllaPSin  S NPLIllaS  S QOTO  IS 

POJTl I l■P5l I I S NPL(2»aS 

CONTINJE 

IPIlMl r.C0.2l  90T0  20 

IPI Pbl 1 I .OE.XLl .OR.PSI 11  .LE.O.OI  OOTO  203 
|FIPb(2l.3E.XL2.0R.P6(2|.LE.0.0l  OOTO  209 
IHIT-IMITal 


00  IS  l-WS 

POJTI I l-PSI 11  6 NPL<2)a6 
CONTINUE 

in  IHIT.E0.21  00  TO  20 
|P|  IHIT.EQ.OI  00  TO  52 

ouestionable  ray 

PRINT  50,H,P 
F0RMATI6E14.Z) 

90  TO  52 
ZS-O.O  S Z9a0.0 
00  21  1-1.9 

ZSaZb«(PIN( IlaPI I 11**2 
£0-ZV*(P0UT( I l-PI 1 11**2 
CONTINUE 

IPIZS.LE.Z91  OOTO  23 
DO  22  lal,9 
ZZ-PINI I I 
PIN!  I laPOUTI n 
POUT! I l-ZZ 
CONTINUE 

ZZaNPLIll  I Z6*Z8 

NPLI  I l■NPLI2l  % Z8aZ9 

NPL(2I*ZP  S Z9aZS 

SLR  ■ S0RTPIZ91  • SORTPIZSl 

DO  51  1*1,3 

Pin  • PIN!  II 

CONTINUE 

RETURN 

END 


SUSRUUTINF  CHECK 

C0MMUN/S3RS/MPI3I ,N|9I ,PI3I ,P0I3I ,BULSI I 3 I , SPOT , I SEE  0 , J SEED, 
• SNT.CST 

COMMON/LATIX/8RIOI4I ,BKPT, I BKP T , SPL T 1 1 • 1 
00  1 I a 1,2 
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A 


r" 


I 


ir(M( I I .NE.O.O)  00  TO  1 
X ■ p( n /aiiiD( n 
tx  ■ xnxr(x> 

iriX.NE.ELOATF( IX) I OU  TO  I 
X • JRANOll ISCED) 

Pill  ■ Pin  * 0.00001 
IP(X.OT.O.S)  Pill  • Pin  - 0.00002 

1 CONTINUE 

IPINI3I .NC.O.OI  90  TO  2 
IFIPI 31 .ST.BKPT)  00  TO  3 
X ■ Pl3)/9RIDI3I 
IX  ■ XPIXPIX) 

IPIX.NC.FI.OATFI  IXI  I 00  TO  2 
90  TO  4 

3 X ■ IPI3)  • BKPT)/ORIOI4l 
IX  ■ XFIXFIX) 

IFIX.NE.FUOATFI IXI ) 00  TO  J 

4 X ■ JR4N31I ISCED) 

PI3)  ■ PIS)  * 0.00001 
1FIX.9T.0.9I  PIS)  ■ PI3)  • 0.00002 

2 CONTINUE 
RCTUHN 
END 

C 

C 

SUBROUTINE  CH099ISURCI 

CONNON/laTCN/JLAT 12. lOUl ,NLATI20) . ICATiNLATT. ISECT. ILATX, ISECTT, 

• PLAriA.lOni.pTARIO.IOU ) .PODV.PBOX. SCPT.NOFL. I H I ST . I PT . N ARC , NPT S 
CONMON/POINT/PSISI.PRI J|,PEI3I,PCISI.0ELR,NFL03 
CONHOn/SORS/NPIS) .HIS) >PISI ,P0I3I«RUL91IS)«SPCT, ISEEO.JSECD, 

• SNT.CST 

CONMON/laTIX/OKI 014) .BAPI, I BXPT , SPUT I 1 9 ) 
C0NNUN/LAT/NX,NT.N2.N3tCT. I SL 
CON  OETCNNINES  if  section  CK0S3IN9 
ISECTT-ISECT*! 

8LRC  ■ 2000000.0 

IFIHI3).SE.O.O.AND.I3ECT.EO.NSECT)OOTO  3 
IFIHIS)  .UE.O.O.ANO.  ISECT.EU.l  ) 00  TO  3 

IFIhIS)  .£0.0.0)  00  TO  3 

IFIhIS) .lT.0.0)  SO  TO  1 
PCIS)  ■ SPlTIISECT) 

90  TO  2 

1 PCIS)  • SPLTIISECT  - 11 
ISECTT  • ISECT  - 1 

2 CONTINUE 

SLRC  ■ IPCI3)  - PSI3) )/HI3l 
PCI2)  ■ PSI2)  ♦ »I2)*9LRC 
PCI  1 ) • PSI  1 ) « HI  1 )*SLRC 

3 RETURN 
END 
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o o 


suanoUTHE  edit 

C0'1iaN/M*NN/MAN(ll0,55»6l,NTE8t 

COHMU^/EVAU/VI(200),N¥l(200),Vl2(200),r*CTHl) 

C0>(MJ  </T'<AC/SLW,XLI,XL2,*L3,NH1TS,DIST(2.1'1I,  t0IST(10|»NCELt 
C0'iMJN/SaRS/wP(3) > M ( 9) >RI 3) ,P0I 3) r BUL9I I 9) > SPCT/ I SEED, J SEED# 

• SST.CST 

100  FORUrdR  I 

101  F0RHAT(Aln,5X,A10,Al,2X,A10,Al,9X,A7l 

102  FOR-iAr(110,A5,E11.4,A2,E11.4,A9,|71 

103  FORMAdSAlOl 
110  FORHATl  Ml  1 

CAfI  • lOH  INCREHENTS  CAP2  ■ 9H*  • *1  CAP3  ■ 2H4* 

CAF4  • FiMO  HAYS  S CAPS  ■ IQH  STAND  DEVS  CAPA  ■ IN 

CAP7  P IDH  VUL  INDES  CAPA  ■ 1 hx  1C  A P »■  1 0H«  > 

CALL  YAH 
1Z*0 

10  4PIT£(6,110) 

NR|TE(6,101)CAP1,CAP7,CAPB,CAPS,CAP6,CAP4 
HP  I TE  I 6,  103)  CAP9,  C AP9,  CAPS),  CAP9, CAPS 

11  17  ■ IZ  A 1 

112  HPITe(S,102)  |Z,CAP2, VI ( 17) ,CAP3,VI2I I7),CAP2,NVI < 17) 

1F( NVI ( 17«1 ) ,EU,0)  90  TO  20 
if(huu(u,50).eq.0)  so  to  10 
SO  TO  11 
20  RETURN 
END 
C 
C 

SUBROUTINE  EVEN(JEXIT) 

COM  BRINSS  PBDY  TO  EVEN  NUMBER  OF  OELR  OR  EXITS  BOX  OR  CROSSES  SECTION 
C0b‘'iin/nann/)'AN(110,55»6),nteST 

CONflUN/LATEN/JLATI2,10U),NLAT(20), I L AT , NL AT T , I SECT , I L A T X , I SEC TT , 

• PLATI4,100),PTAR(6,10U),PBOY,PBOX,SCPT,'JnFL,IHlST,  I P T » N A RC , NPTS 
CORMnN/EVAL/VI(2D0),NVl(20l)),V12(2pn),FiCT(ll) 
C0-)NUN/S9R3/HP(3),M(3),P(3),PU(3),BULSI(3),9PCT,ISECD,JSEED, 

• SNT,CST 

COB^ON/TRAC/SLR,XL1,XL2,XL3,NHITS,DIST(2, 10) ,1DIST«  10),NCELL 
CnBNUS/P0INT/PS(3),PR(3), PE ( 3 ) » PC ( 3 ) , OELR, NFL05 
JEXI r • 0 

7 call  CROSS(SLRC) 

PPEL  ■ PSDY  - OELRaFLOATFINbEL) 

79  ■ DELR  • PDEL 

|F(SLRC.L£.79)  GO  TO  20 

IFISLR  • PiOX.LT. 79  ♦ l.OE-10)  00  TO  2B 

call  LATISIPS, I , J,K) 

call  LOCATEIK, ISECTT) 

1F( hAnI  I , J,X  ) .EQ.U  ) 00  TO  31 
RNl  • URAN31 I ISEED) 
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r 


1 


BEST  mum  cc: 


r>n>/ 
l/i  i 


DO  2 (Z  • t,3 

rE(  IZ)  ■ PS(  IZ)  « Z9*«l(  IZI 
PPIIZ)  • PS(IZ)  * Z9*i4l  |Z)*RN| 

2 CONTINUE 

CALL  UATISIPE, I.J.KI 
CALL  located,  ISECTTt 
IPINANl I, J,KI.EQ.0)  00  TO  3 
call  LATISIPK, I, J,Kt 
call  locateia, isectt) 

IF(MAN( l,j,K).EU.OI  90  TO  3 
call  SCUNE( sc, ISCf I> J,A ) 

SCPT  * SCpT  * Z9*SC*FACT(  I3C)/0ELII 

INK  ■ n 

PPITE(&, 1000)  PS,PE,P,INX 
1000  F0SMAr(9(E12.4,lX),/,|9) 

00  130  IZ  ■ 1,3 
PS( IZ)  ■ PEC IZI 
130  CONTINUE 

PBDY  • 0EL9*FL0ATF(NDEL«1 ) S P80X  ■ PBOX  * Z9 
SOTO  91 

3 call  TBACE(Z9I 
IZ  ■ 0 

9 IZ  ■ IZ  A 1 

IFCDI3T(2, IZl.EO.O.O)  SO  TO  14 
PDEL  • PUEL  ♦ DIST(1,IZ) 

ISC  ■ IDISTC IZ) 

ifcpuel.lt.delr)  00  TO  i; 

UIST(1,IZ)  - DISr(l,lZ)  > POEL  * OELR 
NCELL  ■ IZ 

17  SCPT  • SCPT  A OISTI  1 , I Z)  AD1STI2,  IZIaFACTC  ISO/OCLB 
PBDY  • PBDY  A DIST(1,I4) 

14  PBOX  ■ PdOX  A DISTC I, IZ) 

00  6 IX  ■ 1,3 

PS(IX)  ■ PSCIX)  A OISTI 1, IZ)*N( IXI 
6 CONTINUE 

IFI IZ.LT. NCELL)  90  TO  9 
INX  ■ 12 

NRITtCb, 1000)  PS, PE, N, INX 
IFCPDEl.LT.DELPI  90  TO  7 
90T0  91 

9 CALL  CkOSSISLRC) 

31  IF ( SLNC.lE.ZB)  00  TO  20 

IFISL"  - P30X.LT.Z9  A I.OE'IO)  00  TO  21 
UO  9 IZ  < 1,3 

PEC  IZI  ■ PSC  IZ)  A Z9A«CI  IZI 
9 CONTINUE 

call  LATISCPE,!, J,K) 
call  lOCATEIX, ISECTT) 

IFCNANI I, J, XI. ST. 01  00  TO  3 
DO  62  IZ  > 1,3 
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BEST  AyAILABLE  COPY 


ps(  li)  ■ pe( iz> 

62  CONTINUE 

P90X  • PBOX  ♦ 29 
29  ■ UELR 
INK  • 13 

NHlTfc(6>1000l  PS«PE>N,tNX 
SOTO  8 

20  IPISLH  . P80X.LT.SLRC  t I.OE'IO)  00  TO  2S 
IFLXS  ■ 1 t 2u  ■ SlRC 

90  TO  29 

28  IFLXS  ■ 0 S 29  ■ SLR  « P80X 

29  call  TRACE(Z9) 

INX  • 14 

mR I TLI b. lOUOl  P3,PE,M>INX 
IZ  ■ 0 8 IFLA02  ■ 0 

21  tZ  ■ 12  »1 

IFI UI8T( 2, I 21 .EO.U.OI  SO  TO  24 

ISC  ■ IDISTI IZI 

FUEL  ■ PUEL  ♦ DISTI 1» I2> 

IFIPOEL.LT.DELRI  so  TO  22 
UIST(1,12)  ■ 0IST(1,IZ)  • PDEL  * DELR 
NCELL  ■ 12  S IFLAS2  ■ 1 

22  SCPT  ■ SCPT  ♦ DISTI 1. U)«DlSTt2, IZ)»FACTI ISCI/OELR 
PBDT  • P8DT  ♦ DISTI 1» |2I 

24  PUOX  ■ P80X  * DISTI If  12) 

DO  23  IX  ■ lf3 

PSIIXI  ■ PSIIX)  » DISTI If IZ)*N| IX) 

25  continue 

IFI I2.LT. NCELL)  SO  TO  21 
IFI IFLASZ.EO.l ) SO  TO  bl 
JEXIT  • 1 

IFI  IFLAS.EOfO)  GO  TO  92 

IFINFL33.Ea.3l  NLATT  • 0 

IF  I iLATX.Or . 1 ) SO  TO  4U 

NLATT  ■ NLATT  ♦ 1 S 12  ■ NLATT 

SO  TO  41 

40  12  ■ ILAT 

41  CONTINUE 

NLATIISECTT)  ■ NLATIIStCTT)  ♦ 1 
call  LATENT  I 1 f I2f PC) 

I F I NFi.33.Lr  .3)  00  TO  100 
call  LOADI ISECTT) 

ISECT  ■ ISECTT  1 ILAT  ■ 1 
DO  TU  lOD 

92  29  ■ uELR*FLUATFINUELI 
IFIP8iJT.EU.Z9l  60  TO  50 
KN1*UNAN31 I ISEEO) 

IFI«Nl.9T.PH0y-29IG0  TO  90 
SC  ■ SCPI«0ELR/IP8DT  - 29) 

SO  TO  91 


I 

I 
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so  sc  ■ 0,0  % 00  TO  too 

SI  NOEL  ■ NOEL  * 1 

call  TALurCNUCLfSCflll 
too  SCPT  ■ 0.0 
ketuhn 
ENU 
C 

c 

SU8000T1SE  HCHI 

COHnOr<t/SORS/NP(3l,N<3l>PI3l,EOI3)>BULSI  I3)>SPCT,  ISEED.JSECOf 
• SsTfCsT 

COfIMUM/MI SC/PI >P2> OCR, SCH,NPL01>  NFLS2, sot « $02, 0n1«0m2 
IFINFLOl.EO.alUOTU  6 
CSTHT  ■ JRAN31 ( ISEEDi 
KaR2 

ooro  I 

b xaSORTP ( 3CR**2aOCH**2| 

CSTMT-OCR/R 

1 XI  « 0RAS31 I ISEEOi 

X2  • URAH31 ( ISEEOI 

X|SaXI*«2 
X2SaX2**2 
T*XtS*X2S 
IFIT.OE.t.OIOOTO  1 
CSPMI«( X1S-X2SI/T 
SNPMl"2,U«XI«X2/T 

XI  ■ iJMARJt  ( tSEEDI 

iPIXI.OT.O.S)  UOTO  2 
snphi«>srpri 

2 SRTHT»saRTF( 1 .0-CSTHT»*a» 

API  1 l•SMTHT•SMPH^ 
t<P(2l*SNTHr»CSPMI 
«P(3)«CSTMT 
IFInFL3I.E'J.2)  OOTO  S 
|F(NFL3t,or.3(  GOTO  3 
nP(3I*«RP(3I  S goto  S 

3 IFlNPLGI.OT.S)  OOTO  S 
X 1-A43F ( RP( 2) I 
RPt3l«X|«-»(3l/RP|2l 
IF ( NFL3I .ST .4)  OOTO  4 
«P|2)«Al  % SOTO  S 

4 i.P(2l«.Xl  t SOTO  S 

5 IFINFLSI.ST.FI  OOTO  S 
XIaAPSF I 4P( I ) I 
4P(S>aXl»4PI3)/4P(n 
IFTNFLOI.ST.S)  OOTO  7 
RP||)*AI  1 SOTO  S 

7 4P(l>a«XI 

8 00  10  Ial,3 

P0(|)aR*4P(n  S RP(l)a*NP(|) 
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r 


p( n • Poll) 
to  CONTMUC 

KETUMN 

CNO 

C 

C 

SUORUUTInE  Inx3 

COM  calculates  exit  point  fHOH  CELL 

COMrtON/T^AX/pi ( 31 »P2( 3) ,OOUNO( 6 I . PLl • PL2»PL3. E9 
CONMUN/S3Ri/Np(3),ii|(3)«pl3),POI3),8uLSH3)rSPCT,lSEED,J8CED, 

• SNT»CST 

IPI Ml i ) .CQ.n.O)  so  To  2 

P2(  1 1 • PH 

19  • (PLl  • PI  ( I ) l/WH  I 
P2I2>  ■ P1I2)  * Z9*N(2> 

P2I3)  ■ PI  (31  * Z9*M(3) 

IE(P2( 21 .9T.B0UNU(A)  .0K,P2(2)  .LT.90UND(3) ) 00  TO  2 

|F(P2( 31 .9T,B0UNU(6I .UM.PZl 3) .LT.SOUNOIS) ) 80  TO  2 

SO  TO  » 

2 IFIN(2) .EO.O.O)  SO  TO  3 
P2(2I  ■ PL2 

19  • (PL2  • P1(2I 1/8(21 
P2(  1 1 ■ PI ( 1 ) « Z9*H( 1 ) 

P2(3)  ■ Pl(3)  * Z9*8(31 

IF(P2( 1 1 .ST.B0UNU(2) .0P,P2( 1 I .LT.BOUSDI  11  I SO  TO  3 
IF(P2( 31 .ST.BOUNOlPl .U8.P2I 3) .LT.SOUXDISl ) SO  TO  3 
SO  TO  S 

3 IF(8(3I. EO.O.O)  SO  TO  4 
P2(3I  ■ PL3 

Z9  • (PL3  • PI (3) )/8(3) 

P2( 1 1 • PI ( 1 ) ♦ Z»»8( 1 ) 

P2(2)  ■ PI (2)  * Z9«8(2) 

IF(P2(  t t .ar.UQUHUli)  t)  ,LT,800’iD(l)  ) 00  TO  4 

IF(  P2(  2 1 .ST.SOUND(  4 I .U'<,P2(  2)  .LT.eaUNO(  3)  ) 80  TO  4 
SO  TO  9 

4 CONTINUE 
COM  TRACE  ERROR 

5 CONTINUE 
RETURN 
END 

C 

C 

SUBROUTINE  I80T 

C0MM0N/S0R3/i«P(3) .m(3) >P(3) >P0I3) >8UL3I (3) fSPCT, ISEED.JSEED, 

• SNT.CST 

C0MMUN/Nlsc/Ml,R2>9CR.SCH,NFLSlf NFL02,SDlrS02,QW)>3M2 

XI  ■ JRAN31 ( iseeo) 

)F|R2.3E.R1)  SOTO  1 
C9THT«2.UoXl-l  .0  I SPCT-l.O 
SOTO  2 
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BESrAVAlKJ 


1 X2  ■ SaRrF(R2**2  • Rl**2l/n2 
SPCT  ■ ».0  - X2 
CSTHT  ■ J.O  - Xl»(l,0  • X2I 

2X1  ■ JRXNJH  ISEEO) 

X2  • JRXNJt ( ISEEO) 

X|3«X1**2 

X2S«X2**2 

r-Xt3«X2S 

IF(T.8E.1.0)  00  TO  2 
C3PH1-(X1S-X2S)/T 
SNfiHl«2,U»Xl»X2/T 
XI  ■ jHn'tiH  ISEEOl 

1E(X1.3T.0.5)  GOTO  J 
SRPH1«-3RPMI 
3 T1»«P( n ••2*«P(2»«*2 
T3«l.0«CST-iT»»2 
T4-T3/T1 

IE(T1.3T.1.0E.tO)OOTO  4 
3'iTHT»5aRTF  ( T3I 
«! I l-iRTHT^CSPHI 
2l»3STHT'>3apMI 
«(  3>«CST»1T*'<PI3» 

3070  S 

* T2-SURTF1T4) 

H(l)ar2*(WP(n*)4P(3)*SNPHI4MP(2|*C3PHI>«WP(l)*CSTHT 
«»(  21-T2*  ( >(»(2)  •wPI  3»*SNp»i|«i,p(  1 I •CSPHI  )»WP«2»*CSTHT 
«l 3)»PP( 3) •CSTHT»TI«T2«SNPH1 
5 RETURN 
ENO 
C 
C 

SU9R00TJSE  EXTENT! IL. IP»Pn 
CON  STORES  *ND  RETRIEVES  EXTENT  RATS 
UI NENSIOH  PI  I 3) 

CORNON/nXNN/NXN ( U0,55»4)*NTfST 

C0NnON/LXT£M/JLATI2,l0U),Nu»T(20»,lEXT,NLXTT.ISECT,lEXTX,lSECTT, 

• PlXT(4,100»,PTXR(6, lOU ( ,PbDY, N80X, SCPT, NDEE* IHI ST, IPT, N4RC, NPTS 
CONNON/SORS/XP! 31,x(3»,P(3),POI31 ,BUES1 I 3 ) , SPCT , I 3E E D, J SEED, 

• SNT,CST 

C0nMQNXTRXC/5ER,XE1,XE2,XL3,NHITS,DIST12, 10> , ID1STH0),NCEEE 
IF!  IE.3T.0)  00  TO  2 
P8UT  • PEXI(l,IPI 
PBOX  ■ PLXr(2,IP) 

SER  • PuXII3,IP| 

SCPT  ■ PlXT(4,1P| 

00  I U ■ 1,3 

P! ( IZI  • PTXR( 12, IP» 

• ! 12  > ■ PTXR!  12x3, IPI 
I CONTINUE 

ISCCT  ■ JEXT<l,IPt 
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r f ■ r-.  ; p 

rAiuiULL 


NOEL  ■ JLAri2,IP| 

30  TO  to 

2 PL*T( I, IP)  • PBOr 
PL*Tt2,IP)  ■ PHOX 
PLATI IP)  ■ SLP 
PL*T( 4, IP)  ■ 3CPT 
DO  J U ■ 1>3 

PTAPI  1 1.  IP  ) ■ PHIZ) 

PTAPl IZ  + 3, IP)  ■ Ml  IZ) 

3 CONTINJE 
JUATtl.IP)  ■ ISECTT 
JLAT<  2, 1 P)  m nDEL 

10  PErONN 
END 
C 
C 

SUaPOUIlNE  LATISIPI , I , J,K) 

DMENSIdn  pI(3)»pH(3) 

C0N«ON/saPS/-P(3) ,«((3) f P(3) ,P0I3) ,9ULSI (3) »8PCT, ISCED, J SEED 
• SNT.CSr 

CONHON/L ATI X/GPI U( 4 I,8APT,I9XPT,3PLT(19) 
CONnOM/LAT/NXjNY.NZ.NSfcCTjISL 
CONPUN/ I vOfX/ 1 1 , J1 ,K1 
CON  locates  cell  I,J,A  POP  POINT  P 
DO  13  1Z«1,3 
PHI  U)«PI  ( IZ) 

13  CONTINUE 
IP  IZ-O 

12  IZPlZPl 

IP|  U.LT.3)  30  TO  102 
IFIPI ( 3) .ST.aKPT)  00  TO  103 

102  xml’ll  IZ)/GHIU(12) 

30  TU  104 

103  x3<PT»PL0ArF( IBKPT) 
xaXBKKr*|PI(3l-dXPT)/3H|0(4| 

104  XaXPlXE(X) 

IFIX.LT.U.O)  K-K>1 
KK'K 

IFIHI IZ) .or.o.o)  KK-KPl 
X<»FLOATF I XK ) 

XXX»X-XK 

IFIPI IZ).3T.0.0)  XXK--XXK 
IF(XXK,ai.l.0E-12l  80  TO  107 
IFI  U.LT  .3)  SO  TO  100 
IF ( KA  .LT  . I 3KPT  ) SO  TO  lOS 
PI(3l»aKPT«(XK-X8KPT)«8pID(4) 

PHI 3)"P| I 3)  *80  TO  3 
ins  PI  I U)aXX«3PlDI IZ) 

PHI IZ)«PI I IZ)  1 80  TO  3 
SOTO  ll,2,10),IZ 
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1 !■!(  ISO  ro  12 

2 J»K  ISO  ro  12 


0 CQnTINjE 

00  4 I Z • 1,3 

i>niz>  • i»i(izi  ♦ o.oooooi»i4(  izi 

4 CO'^TINJE 
00  TO  19 

10  00  U I Z ■ 1 »3 
f’U  UI*P-i(  III 

11  CONTINUE 

I»ltllJ«J*l»KnK»l 
ll-l  >Jt>JSKl«K 
IF( i.yi.Hxi  I ■ I . J 
IF(  I.LT.D  1 • I ♦ 1 
iFiJ.or.xri  J ■ j - 1 
irij.Lr.i)  J • j * I 
lF(<.Or.'(ZI  4 ■ K . I 
If  14.V.T.I)  K ■ K ♦ I 
HETJKN 
ENU 
C 
C 

SOBROUtl^E  LOIDtISECT) 

C L04US  UtrH»LlTr  UESCWIf'TION  INTO  CORE 

C ISECT  EyJ4t.5  section  nuhSER  TO  SE  L 0 A DE  D I RECOR  0 I 1«3 

C NSiCI  tOJ*LS  NUNStN  OF  3ECII0NS  (31 

C NNEH  ejj*l.s  NilHSEN  OF  vOlU  SLICES  IN  L*ST  SECTION  (Tl 

C01nv)N/«»nn/«I4N(110,55>6),NTE3T 
C0il0N/L*T/NL,NY,NZ,NStCT,l3L 
DlAENalON  s»he(4) 

c oefine  uisc  files  to  be  utilized 

OITA  0|3(;2/|aHS5024  KNAN/ 

OATA  F1LE2/10HBODT  PARTS/ 

C OPEN  TNE  DISC 

call  UISCDLI J,NAME,DISC21 
IFINTEST.EJ.U  go  TO  999 
call  0I3CR0(FILE2, ISECT, I, NANI 
30  TO  1000 

909  call  TEST(ISECT) 

1000  RETURN 
ENU 
C 
C 

subroutine  LOCATEIKP, IS) 

C0NNUn/LATII</SRI1)(4),84PT,  I 8<P  T , SPL  T T 1 9 1 
CONNJN/LATYN*,NY,NZ,NStCT,ISL 
CON  calculates  nan  K FRON  POINT  K 
IS«0 

21  |3«IS«I 

IFIKP.OT, IS'lSLIOOTO  21 
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BKI  COPY 


■« 


KP  ■ KP  . ISL«< IS  - n 
PETUPN 
END 
C 
C 

subroutine  SCqREISC, I3U, I , J,kI 
CONNON/tUNV/NANf  I iu>ss>6i,iNrEsr 
ISC  ■ NAN(1,J,K) 
sc  ■ EtOAT"!  ISCI 

NETURN 

END 

SUBROUTINE  SORC 

CONNUN/latEN/ JLAT(2,  10U),nLAT(20>,IUAT»NLATT»I8ECT,|LAIX,ISECTT# 

• PLAT(A,lon),pTAR(6,10U|,pHOY,PBOX,SCPT,NDEL>IH|ST,IPT>NARC,NPT3 
C0iN0N/3  3R9/>iP(3),>»(3).P(3),P0(3l,BULSI(3),  SPCT,  I8EE0»  JSEED, 

• SNT.CST 

C0iN0'</TRAC/3LR»»Ll,XL2,XL3,NHITS,DlST(2,10l,  I 0 1 S T ( I 0 I » NCEUL 
C0N''(iN/MISC/Kl>R2>UCR.UCH,Nf'LSl>NELS2>SDl>SU2,SWlfyw2 

iF(  ihist.gt.1)  go  to  U 

IE( NFL31 .GT. 1 I GO  TO  8 

ifinpts.gt.i ) go  to  2 
UO  1 1 ■ 1,3 
P( I)  ■ PO( I ) 

1 continue 

30  to  4 

2 DO  3 I ■ 1,3 

P(  I ) ■ PIARI I , IPTI 
<t  ( I ) ■ PTAR  ( I >3,  IHT  > 

3 CONTINUE 

; 4 CONTINUE 

IF1NKL32,0T,2)  00  TO  13 
IF(NFL32.GT.n  SO  TO  5 
R2  ■ (<'(  I ( • 0.S*XLU«»2 

R2  ■ R2  ♦ (P(2)  - 0.5«*L2)««2  i 

R2  ■ R2  ♦ (P(3)  - 0.5«AU3I«*2  '! 

R2  • SJRTF(R2»  V 

30  TO  11  { 

5 CONTINUE  I 

RP  ■ 0,0  i 

0 0 6 I • 1 , 3 y 

R2  • K2  ♦ (P(II  - BULSl(m*«2 

6 CONTINUE 

IN2«SUHrF(R2)  |i 

UO  7 I ■ 1,3  f| 

xPI  I ) • (BJLSI  III  - Pin  I/R2  ! 

7  CONTINUE  jl 

30  TO  12  I 

8  IF|NFL31,9T.8I  GO  TO  9 ( 

I call  NENI  { 

[ GO  TO  10  } 


rn?Y 

I si  iki  L/ 


c 

c 


c 

c 


c 

c 


3U8RoUTI«IE  TALLr(  IP,SC,  ISC) 

COn«0-</H*>*'J/nANI  U0,55»6)  / NTEST 
C0'iNUi'</E»*U''Vlt20U),NvJ  1200) , vl2l200),F*CTt  III 
SC  ■ SC«EACTIISC) 

VI ( IP)  ■ VI ( IP)  ♦ SC 
VI2(IP)  ■ VI2( IP)  * SC*SC 
i NVl  I IP)  ■ WHIP)  ♦ 1 
RETURN 
END 


SUBROUTINE  TESHISeCT) 
C0NMUN/M*NV/«AN( U0«  3Sr 6) ,NTEST 
00  12  I ■ Hlio 

UO  12  J ■ },5S 

00  12  < ■ 1,6 

HAN( |,J,<)  ■ 0 

12  continue 

DO  1 1 -SO, 60 

DO  I J -IjSS 

00  1 K -1,6 

NAN( I, J,«)  ■ 1 - SO 

1 CONTINUE 

00  2 J >24,32 

00  2 I >1,110 

DO  2 K >1,6 

nAN( I , J, < ) > J - 23 

2 CONTINUE 

IF( IbECr.Sr.8)  00  TO  5 
IF( IStCT.LT.F)  GO  TO  5 
DO  3 1 > 1,110 
00  3 J > 1,55 
00  3 6 > 1 ,6 
nANI 1 , j, < I > 10 

3 continue 

5 RETURN 

end 
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con 


con 


SUBROUrME  TK*CE(Rl 

THACivs  ThbduQh  cells  Ano  STORES  SCORES  AND  UISTANCCS 
C0nnuN/bJRS/*PIJ)>"(3)>i>(3l,P0<3lfBJLSI(3lf  SPCT« ISEED* JSCCD, 
* SNT,tST 

C0nMUN/T»*x/Pl(31.P2(3l ,HOU>JD( 6 I , PL I « PL2*  PL3» 19 
C0nHjti/TRiC/SLR.XLl,XL2,xL3,NHITS,DIST(2«t0J, | U 1 5 T ( 1 0 1 > MCCLL 
C0n10n/L»TIX/GRID(4(,B^PT,lBKPI,SPLT^l9) 

COirtjn/POI Ml/ps I 3 ) » PR ( 0 ) ,pt ( 3) »PC ( 3 » » BtLR* RElOS 

COnnON/l sDEX/1 I, Jl,Kl 

initialue  arrays 

00  11  J ■ 1,9 
|OIST( Jl  ■ 0 
00  11  I ■ 1,2 
PI3T< 1,J)  ■ 0.0 

11  CONTINUE 

00  12  J ■ 1,3 
PKJI  ■ PSIJI 

12  CONTINUE 

NCELL  • 0 S PCELL  ■ 0.0 
10  NCELL  • NCELL  ♦ 1 
call  L*TIS(P1,1,J,KI 
PLl  ■ G«lD(l)«FLOATF(n) 

PL2  ■ 3R10(21*FL0ATFIJI| 

IFdC.GI.  IR<PT|  go  to  1 
PL3  ■ 3R10(3)«FL0ATF(KI» 
im  ■ 3RID(3I  S go  to  2 

1 PL3  • BXPT  ♦ QR1D(4(*FL0ATF(K1  - IBKPTl 
ZNZ  ■ SR1DI4I 

2 CONTINUE 

B0UN0(2I  ■ PLl 

B0JNU(4)  ■ PL2 

aOJNOlS)  ■ PLJ 

aOJNOl II  ■ PLl  • GRtDI  1 I 

aOUNUli)  ■ PL2  - GRIDI'Zl 

BOUNUIS)  ■ PL3  - ZNZ 

IF  I «(  1 I .LT.O.i)  I PLl  ■ PLl  • QRID(l) 

IFIxlZI .uT.O.OI  Pl2  ■ PL2  - GR1D(2) 

IF(«( J) .lT.0.0)  PL3  ■ Pl3  • ZNZ 
call  1NX3 

call  LOCATEIR, ISECTTI 
call  SCOREISC, ISC, I , J,R) 

OtSTIZ, NCELL)  ■ SC 

IDISTINCELLI  ■ ISC 

IFIPCELL  * Z9.GE.K)  3U  TO  14 

01ST(1,NCELLI  • Z9 

PCELL  • PCELL  ♦ Z9 

00  13  IZ  ■ 1,3 

PlI  IZI  - P2MZ1 

13  continue 

IFIR  - PCELL. LT,10E-12I  GO  TO  7 
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90  TO  10 

14  OISTd.NCELL)  • R • PCCLL 
7 COnTMuE 
WETUHM 
END 

c 

C 

SUaROUTHE  TR*CK«H.I 
C0aN0N/>i*NN/M*N(ll0,5S»6»,NTE8T 

Cn>iMjn/L*TEN/ JuAT  12,  lOUl  ,NI.*T«20>  , ILtTf  NL4TT.  ISECT.  1L4YK,  I9ECTT* 

• PLAri4,|OOI,PTAH(6,10U),,>dDY,HaOX,SCPT,NUEL,IH|ST,IPT,NARC,NpT3 
CONMON/E 1(2001, NVi(20a|,vi2(200l,PACT(in 
C0NndN/33RS/»P(3),H(3(,P(3l,P0|3), BUUSI ( 3) , SPCT, 1SEED,JSEED, 

• SNT.CST 

C0oihJN/T4AC/SLR,XL1#XLZ,XL3,NMITS,P1ST(2, 10), lOISTl 10|,NCELU 
C0’1Ndi«(/P3IsT/PS(J),PN(3»,Pe(3),  PC  1 3 ) , DELS,  NFL03 
COinON/L*T  I X/ONll/(  4) , aAPT,  I BKPY  , SPlTl  19) 

COinoN/lSDEX/Il,Jl,Kl 

COnHUiM/LDOE/NLOAO 

cnti  TRACKS  IxiriAU  «AT8,LAIENTS,0R  RAYS  returned  prop  even 
SU  JEX I T ■ D 

IPlIL.EQ.l)  00  TO  120 
call  LATENTIO, ILAT,P3) 

■R1TEI6,523)  ILAT,NLATT,IseCT,ISECTT,JEXIT,NDEL 
525  PORMATl 10(2X, 15) ) 

JUATl  1, ILAT)  ■ 0 
NLATIISECT)  ■ NLAT(ISECT)  • 1 
IFIPUOT.EO.0.0)  80  TO  2 
call  EVENIJEXIT)  S 00  to  4 

120  DO  I U ■ 1,3 
PSl IZ)  • PI IZ) 

1 CONTINUE 
INX  > 5 

NRI TE 1 6, 1000)  PS, PE, N, INX 
1000  PONNAT I 91 E12.4, IX ) , /, IP) 

PBDY  ■ O.O  S POOX  ■ 0.0  $ NDEL  ■ 0 

2 continue 

4 IF! JEXIT.Ea.l)  00  TO  70 

132  call  l.ATlStPS,l,J,K) 
call  L3CaTE(K, ISECTT) 

IFI IseCT.EO.lSECTT)  UO  TO  43 
CON  STORE  LATENT 

UO  25  IZ  ■ 1,3 
PCI  IZ)  ■ PSl IZ) 

25  continue 

133  |F(NFL33.Ea.3)  NLATT  ■ 0 
in ILATX.GT.I I 00  TO  40 

NLATT  ■ NLATT  ♦ 1 f IZ  ■ NLATT 
80  TO  41 
40  IZ  - ilat 
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41  CONTINUE 

NL4T(IseCTTI  ■ NLATCIStCTTI  ♦ I 

CALC  UATC>jr(  1,  U.PO 

IF »NFL33.LT,3I  00  TO  70 

call  UQAu( ISECTT  I 

ISECT  - ISECIT  I ILAI  ■ I 

00  TO  SO 

CON  ENo  UF  LATENT  STONAOE 
43  CONTtNoe 

call  C40SSISLRC) 

1F(SLHC.LE.DElR»  80  TO  20 

IF(sln  • P^OX.lT.UELR  * I.OE-IO)  00  TO  30 

1 F t riAN(  I , J,  K I .Eu.o  » 00  TO  7 

RNl  ■ JRAN31 ( I SEED) 

DO  S2  17  ■ 1,3 

PE( Ul  « PS( IZI  ♦ DELRoW( 17) 

PR)U>  ■ PS<I71  ♦ OELR»i»(  I71«RN1 
S2  CONTINUE 

call  LaTISIPE, 1, J,K) 
call  L0CATC(4, ioectti 
IF) iAN( I , J,K ) .EO.O)  GO  TO  3 
CALL  L4TIS(PR, I, J,KI 
call  L0C*TE)<, ISECTT) 

IF  ( N* ')(  I , J,K  1 .EO.O  ) 00  TO  3 

call  SCORE! SC, tSC, 1 ,3, A) 

NOEL  N NOEL  ♦ I 

PSUr  • 0£LR»FLOATF(NDEL) 

P90A  ■ PiOX  ♦ DELR 
call  T4LlY(NUEL,$C, ISC) 

IN*  ■ 1 

NRITtlS.lDOO)  PS,PE,M,INX 
00  130  17  • 1,3 
PS(  171  ■ PE( 17) 

130  CONTINUE 
30  TO  43 
3 call  TRAJEIDELR) 

00  3 17  • I.nCELL 

IF! 01 JT! 2, 171 .EO. 0.0  I 00  TO  14 

|SC  • IDlsr(l7) 

SCPT  • SCPT  ♦ DI  SI  ( 1,  I 7 l•DIST(^,  I7I»FACT(  MO/OEUR 

paur  ■ Psoi  ♦ DISK  1,  17) 

14  pao*  • pap*  ♦ UIST! 1, 17) 

5 CONTINUE 
|Na  - 2 

nRI  T£(  b,  moo  I PS, PE, 4,  IN* 

00  6 17  • 1,3 
PS! 17)  • PEI  17) 

6 CONTINUE 

CALL  EVENIJEXIT) 

00  TO  4 
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S CAUL  C40SSISLRC) 

1 iR(SLRC.Le.DELR)  SO  TO  20 

IFtSUX  • RSnx.LT.DELM  * I.OE'IU)  00  TO  30 
DO  9 IE  ■ 1,3 

RE(  Ul  ■ PS(  12)  * OELR*M(  |Z) 

9 COMTi'XUE 

C«LL  LxriSIPE, I.J,K) 

CALt  L0CATE(X, ISCCTT) 

IP ( MANt I > J, K > .OT.O)  00  TO  3 

61  DO  62  12  ■ 1,3 
PS( 12)  • PEI  12) 

62  CQxTi'lJE 

PBOX  ■ P80X  ♦ DEL* 

00  TO  a 

20  IPISlM  • PSOX.LT.SLRC  * I.OE'IO)  00  TO  30 

call  traceislmci 
DO  21  12  ■ 1,NCELL 
IF  1 DIST I 2, I 2) .EO.U.O)  SO  TO  24 
ISC  ■ lOISTI 12) 

SCPT  ■ SCPT  ♦ DISTI 1 , 1 2)*01ST(2, 1Z>*PACTI |SC)/DELR 
PBUY  ■ PSDT  ♦ DI$r(l,U| 

24  PSOX  ■ PSOX  * UlSTI 1, 12) 

INX  ■ 3 

kRITE(6, 1000)  PS, PE, H, INX 

21  CONTINUE 
SO  TO  133 

30  call  TRACEISLR  - PSOX) 

00  31  12  ■ l,NCELt 
tF(DISTI2, 12).EO.O,0)  SO  TO  34 
ISC  ■ lOlSTI 12) 

SCPT  ■ SCPT  ♦ DlSTIl,  U)«DlSTia,  t2)«PACTnSC)/0ELR 
poor  ■ Psnr  ♦ oiSTi i, 121 
34  PSOX  • PSOX  ♦ DISTI  1,  121 

31  continue 

INX  ■ 4 

>iRlTEl6,tnOO)  PS, PE, N, INX 
28  ■ UELM^FLOATFINDEL) 

IPIPS0T.E(1,28I  60  TO  70 
MN1-UNAN31 I ISEEDI 
IP|RNl.3T,P8UY-28)80  TO  70 
SC  ■ SCPr»OELR/l P8DY  - 28) 

51  NDEL  • NOEL  ♦ 1 

call  TallYINUEL,SC,11) 

70  SCPT  • 0,0 
100  CONTINUE 
METUHN 
END 
C 

c 

SUSROOriNE  TRANS 
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CnM  T«»*»JSL*TES  POINT 

CONHlj't/If<AC/SLR.XH.XL*,XL3,NHITS,DI»T(2,lO).IDI8T«IO»*MetLL 
CON«UN/SORS/i»P(3I ,N(3» ,P(3) ,pnt3) , BUU6I ( 3 > . SPOT , I SCEO* 3 SEED, 

• SNTjCSr 

C0nmUN/NIsC/R1,R2.0CR,UCM,NFL81,NPL02,SD1,SD*»9W»*#N* 
p(  1 )»p(  1 i^a.i^xn 
P(2)»P(2)*0,6»XL2 
IFlNKLSl.EJ.bl  GO  TO  too 
P(3)"P(3U0,9«X1.3 
too  RETUNN 

END  4 

C 

C j 

SUBROUTINE  V*R 

COnMUN/mann/maNI I 10,5S>6I f NTEST 
COIinon/E»al/VI(2QU>,Nv1(200),V12(200).P*CT(1II 
IZ-0 

11  IZ-IZ*I 

IF(NVl ( IZI.EU.O)  GOTO  20 
IFiNV I ( I Z I .3T. 1 I GOTO  12 

VI2(  1Z)-R.S)9  S GOTO  It  ; 

12  l■NVI(IZI 
Z8»FL0ATF(  I » 

Z9>VI  ( IZ)/Z8 
VI(  IZIaZ9 
29-Z9*Z4*2B 
Z9>V12( tZI<Z9 

ir(Z4.LT.O,O.AND.Z9.QT.<1.0E-06)Z9>0.0 
|F( Z4l 21, IB, 16 
21  NRtTEI6,100IZ4 

too  F0RMAT(3UX, 'THE  VALUE  OF  Z9  I 3 • , 2X i E 1 2 . 9 1 
30  TO  11 

16  Z9-Z9/(ZS«1.0) 

Z9«3NRTF(E9I 
V12( 1Z)*Z9 
GO  TO  11 
20  RETURN 
ENO 
C 
C 

• COMPILE  UISC,NRAN31,ALL 

• CORPILE  U1SC,URAN31,ALL 

• LIST 

• DATA 

BULSl  CHECX  8IVAR  DISTRIB 
11111111  11113111 

1 I I 


r 


DISTRIBUTION  LIST 


No.  of 
Copies 


Organl zation 


No . of 
Copies 


Organ i zation 


12  Commander 

Defense  Documentation  Center 
ATTN:  DDC-TCA 
Cameron  Station 
Alexandria,  VA  22314 

1 Director  of  Defense  Research 
and  Engineering  (OSD) 
Washington,  DC  20305 

1 Director 

Weapons  Systems  Evaluation  Gj) 
Washington,  DC  20305 

1 Chairman 

Defense  Science  Board 
Washington,  DC  20301 

1 Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN;  DRCDMA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 Commander 

US  Army  Aviation  Rescarcli 
and  Development  Command 
ATTN;  DRSAV-E 
12th  and  Spruce  Streets 
St.  Louis,  MO  63166 

1 Commander 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1 Commander 

US  Army  Electronics  Command 

ATTN:  DRSEL-RD 

Fort  Monmouth,  N.J  07703 


1 Commander 

US  Army  Missile  Research 
and  Development  Command 
ATTN:  DRDMl-R 
Redstone  Arsenal , AL  35809 

1 Commander 

US  Army  Tank  Automotive 
Research  fi  Development  Cmd 
ATTN;  DRDTA-RWL 
Warren,  MI  48090 

1 Commander 

US  Army  Mobility  Equipment 
Research  Development  Cmd 
ATTN:  Tech  Docu  Cen,  Bldg.  315 
DRSME-RZT 

Fort  Belvoir,  VA  22060 
1 Commander 

US  Army  Armament  Materiel 
Readiness  Command 
Rock  Island,  IL  61201 

5 Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN;  DRDAR-LCS,  Mr.  G.  Gaydos 
DRDAR-S^S,  Mr.  Mancini 
DRDAR-SCN,  Mr.  Malinoski 
DRDAR-SCH,  Mr.  H.  Kahn 
DRDAR-SCA,  Mr.  Bornheim 
Dover,  N.l  07801 

1 Commander 

US  Army  Harry  Diamond  Labs 
ATTN:  DRXDO-Tl 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 
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Copies  Organization  Copies  Organization 


2 Conunander 

US  Army  Materials  G Mechanics 
Research  Center 
ATTN:  Dr.  R.  Frost 
Mr.  E.  DeLuca 
Watertown,  MA  02172 

2 Commander 

US  Army  Natick  Research  and 
Development  Center 
ATTN:  Mr.  T.  Keville 
Mr.  M.  Donomee 
Natick,  MA  01760 

] Commander 

White  Sands  Missile  Range 
ATTN:  STEW-TE-PS,  Mr.  Findley 
White  Sands,  NM  88002 

1 Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL  (Tech  Lib) 

Wliite  Sands  Missile  Range 
NM  88002 

1 Commander 

US  Army  Training  and  Doctrine 
Command 

ATTN:  MAJ  R.  Young 
Fort  Monroe,  VA  23651 

i Commander 

US  Army  Infantry  School 
ATTN:  MAJ  B.  Davis 
Fort  Banning,  GA  31905 

1 Commander 

US  Army  Armor  School 
ATTN : Armor  Agency- 
Fort  Knox,  KY  40121 

1 Commander 

US  Army  Field  Artillery  School 
ATTN:  Field  Artillery  Agency 
Fort  Sill.  OK  73405 


1 Commander 

US  Army  Aviation  School 
ATTN ; Aviation  Agency 
Fort  Rucker,  AL  36362 

1 Infantry  Research  & Dev  Ln  Ofc 
ATTN:  COL  John  0.  Ensor 
Fort  Benning,  GA  21905 

1 The  Surgeon  General 
Department  of  the  Army 
Washington,  DC  20310 

1 HQDA  (DAMO-CSS,  MAJ  Rogers) 
Washington,  DC  20310 

1 Chief  of  Naval  Operations 
ATTN:  OP-03EG 
Department  of  the  Navy 
Washington,  DC  20350 

1 Commander 

US  Naval  Air  Systems  Command 
ATTN:  AIR-604 
Washington,  DC  20360 

2 Commander 

US  Naval  Weapons  Center 
ATTN:  Code  753 
Code  12 

China  Lake,  CA  93555 

1 Commander 

US  Naval  Medical  Field 
Research  Laboratories 
Camp  LeJuene,  NC  28542 

1 Commander 

US  Naval  Research  Laboratories 
ATTN:  Code  8434,  Mr.  Ferguson 
Washington,  DC  20375 

1 Conmandant 

US  Marine  Corps 
Washington,  DC  20380 


L. 
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Copies  Organization  Aberdeen  Proving  Ground 


1 

Director 

Cdr,  USAEA,  ATTN: 

MCDEC 

DRDAR-CLB-B,  Mr.  L.  Sturdivan 

ATTN:  C,  Gnd 

Cbt  Div 

Dir,  USAMSAA 

Quant i CO,  VA 

22154 

ATTN:  Dr.  Sperrazza 

DRXSY-GI,  Mr.  E.  Heiss 

1 

Director 

DRXSY-M,  Mr.  R.  Pollard 

MCDEC 

Dir,  USAUEL 

ATTN:  MAJ  D. 

Crowe 

ATTN:  Mr.  B.  Corona 

Quant i CO,  VA 

22134 

Marine  Corjis  Ln  Ofc 

1 

AFATL  (DLRW) 

Eglin  AFB,  FL 

32542 

1 

US  Department 

of  Treasury 

Training  Center 
ATTN:  ^^r.  IV.  McClarin 
Powdermill  Road 
beltsville,  MD  20705 

1 Federal  Bureau  of  Investigation 
ATTN;  Mr.  C.  Cunningham,  Rm  7410 
W.ishington,  DC  30405 

1 International  Association  of 
Chiefs  of  Police 
Caithcrsburg,  MD  20760 

1 Rockwell  International 

Los  Angeles  Aircraft  Division 
ATTN:  Mr.  K.  Dotseth 
International  Airport 
Los  Angeles,  CA  90009 
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